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INTRODUCTION 

The  ultimate  objective  of  the  Northrop/AFWAL  program  in  Reference  1-1 
is  to  develop  an  analytical  methodology  to  predict  the  strength  and  lifetime 
of  bolted  composite  and  metallic  structural  components.  The  development  of 
a  strength  analysis  for  plates  bolted  by  many  fasteners  involves:  (’.)  an 
analysis  that  computes  the  distribution  of  an  applied  load  among  various 
fasteners  in  a  bolted  plate;  and  (2)  a  strength  analysis  that  can  be  applied 
at  every  fastener  location.  The  second  analysis  addresses  a  segment  of  a 
bolted  laminate  that  includes  a  single  fastener,  and  is  the  subject  of  this 
report.  Figure  1-1  presents  a  schematic  representation  of  all  the  analytical 
steps  involved  in  the  strength  prediction  of  bolted  structural  components  . 

At  the  initiation  of  the  Northrop/AFWAL  program  (Reference  1-1),  BJSFW 
was  the  only  available  computer  code  that  performed  a  strength  analysis  of 
a  bolted  laminate  that  transferred  loads  via  a  single  fastener.  Reference  1-2 
presents  the  details  of  the  analysis  that  is  performed  by  the  BJSFM  code. 

The  BJSFM  analysis  can  be  used  to  approximately  predict  the  strength  of  bolted 
joints  when  (1)  the  fastener  loads  have  already  been  computed  using  other 
analyses,  (2)  the  fasteners  are  not  too  close  to  one  another  or  to  a  neighboring 
cut-out,  (3)  the  fasteners  are  rigid  and  their  displacement  with  respect  to 
the  plate  is  uniform  in  the  plate  thickness  dlreclton,  and  (4)  the  first  ply 
failure  at  a  fastener  location  essentially  precipitates  joint  failure. 

The  conditions  mentioned  above  severely  restrict  the  application  of  the 
BJSFM  code  as  an  analytical  tool  in  the  design  of  bolted  structures.  An  im¬ 
proved  strength  analysis  was  developed  in  the  referenced  Northrop/AFWAL  pro¬ 
gram  to  overcome  most  of  the  limitations  in  the  BJSFM  analysis.  Details  of 
this  Improved  strength  analysis,  programmed  to  be  the  SASCJ  (Strength  Analy>:;i.8 
of  Single  Fastener  Composite  Joints)  computer  code,  are  presented  in  the 
following  sections. 
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From  an  analytical  standpoint,  the  major  limitations  of  the  BJSFM 
analysis  are: 

(])  It  assumes  the  planform  dimensions  of  the  bolted  laminate  to 
be  Infinite  (L  and  W  -*<”  In  Figure  1-2), 

(2)  It  circumvents  the  bolt/lamlnate  contact  problem  by  assuming  a 
cosinusoidal  bearing  load  distribution  over  half  the  hole  boundary.  It  also 
ignores  the  effect  of  friction  between  the  fastener  and  the  laminate  (see 
Figure  2-1). 

(3)  It  does  not  account  for  the  effect  of  many  joint  parameters 
that  render  the  two  dimensional  analysis  inaccurate.  Load  transfer  in  a 
single  shear  situation,  for  example,  is  assumed  to  be  equivalent  to  a 
double  shear  situation.  The  fastener  is  also  assumed  to  be  rigid  in  shear 
and  bending,  and  the  effects  of  fastener  torque  are  ignored. 

(4)  Joint  failure  is  assumed  to  be  precipitated  by  the  first  ply 
failure,  each  ply  is  assutued  to  be  linear  elastic  to  failure,  and  ply 
failure  modes  are  not  predicted.  Consequently,  a  comparison  between  pre¬ 
dicted  joint  strengths  and  measured  values  exhibits  a  poor  correlation  in 
many  situations, 

The  SASCJ  code  discussed  in  this  report  overcomes  all  the  above 
BJSFM  limitations  with  the  exception  of  item  (2),  In  overcoming  these 
llmitaticns,  an  Improved  two-dimensional  plate  analysis,  a  fastener  analysis 
and  a  progressive  failure  analysis  were  developed, 

A  brief  summary  of  the  program  contributions  to  the  strength  analysis 
of  bolted  composite  and  metallic  plates,  via  the  SASCJ  code,  is  presented  in 
Table  1-1,  These  accomplishments  represent  a  significant  improvement  over 
the  state  of  the  art  at  the  program  Initiation  stage  (the  BJSFM  computer 
code),  and  aided  the  development  of  a  validated  strength  analysis  computer 
code  (SASCJ), 

Section  2  describes  the  two-dimensional  analysis  that  accounts  for 
the  Influence  of  finite  planform  dimensions  on  the  solutions  for  a  plate 

bolted  by  a  single  fastener.  Section  3  presents  the  fastener  analysis  that 
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Table  1-1.  SUMMARY  OF  PROGRAM  CONTRIBUTIONS  TO  THE  STRENGTH 
ANALYSIS  OF  PLATES  BOLTED  BY  A  SINGLE  FASTENER 


Prior  to  Program  Initiation 

Program  Accomplishments 

Available  Strength 
Analysis  Computer 
Code 

BJSFM 

SASCJ 

Bolted  Plate 
Geometry 

Planform  dimensions  are 
assumed  to  be  infinite. 

Actual  (finite)  planform 
dimensions  are  accounted 
for. 

Bolted  Laminate 
Layup 

Analysis  is  independent  of 
(ignores)  laminate  stacking 
sequence. 

Analysis  accounts  for 
the  actual  laminate 
stacking  sequence. 

Fastener  Properties 

Fastener  is  assumed  to  be  a 
rigid,  frictionless  pin. 

Its  bending  and  shear  effects 
are  ignored. 

Fastener  flexibility 
(bending  and  shear) 
effects  are  accounted 
for. 

Fastener  Torque 

The  effect  of  fastener 
torque  is  unaccounted  for. 

Fastener  torque  effects 
are  included  in  the 
fastener  analysis. 

Load  Eccentricity 

Analysis  only  represents 
a  double  shear  load  trans¬ 
fer  situation 

Analysis  differentiates 
between  single  and 
double  shear  situations. 

Strength  Prediction 

Prediction  is  based  on  first 
ply  failure,  and  demonstrates 
a  poor  correlation  vlth 
measured  strengths. 

Prediction  is  based  on  a 
progressive  failure 
procedure,  assuming  two- 
stage  ply  failures,  and 
agrees  well  with 
measured  strengths. 

accounts  for  fastener  flexibility,  load  eccentricity  and  fastener  torque. 
Section  A  describes  the  developed  failure  procedure  that  assumes  nonlinear 
(bilinear)  ply  behavior  and  predicts  the  progression  of  local  (ply)  failures 
until  the  joint  cannot  carry  any  additional  load  (joint  failure).  Section  5 
demonstrates  an  excellent  agreement  between  SASCJ  predictions  and  test  results 
from  Reference  1-1,  validating  the  improved  strength  analysis. 


SECTION  2 


TWO-DIMENSIONAL  ANALYSIS  OF  A  FINITE  LAMINATED 
PLATE  BOLTED  BY  A  SINGLE  FASTENER 

As  discussed  in  Section  1,  a  two-dinensional  analysis  of  a  finite 
'^■olted  laminate  is  a  primary  requirement  in  the  development  of  an  analysis 
".hat  predicts  the  strength  of  bolted  laminates.  In  the  following  sub-sections, 
a  brief  description  of  the  two-dimensional  analysis  developed  in  Reference  2-1 
is  presented  along  with  sample  predictions. 

2. 1  Governing  Equation 

The  two-dimensional  stress  field  in  a  finite  bolted  plate  is  expressed 
in  terms  of  the  Airy  stress  function  F(x,y)  that  automatically  satisfies 
equilibrium  equations  everywhere  in  the  plate  domain.  The  corresponding 
displacement  solution  satisfies  compatibility  requirements  when  the  following 
equation  is  satisfied  by  the  stress  function: 

®22^’xxxx  “  ^^26^*xxxy  ^^^12  ®66^^’xxyy  ~  ^^16^*xyyy 


a,  ,F, 

11  ’yyyy 


(2-1) 


where  a^^  are  laminate  compliances  as  defined  In  Reference  2-2,  Equation  2-1 
is  the  governing  equation  for  the  problem  of  interest, 

2.2  Complex  Variables  Approach 

A  complex  variables  app’^oach,  described  in  Reference  2-3,  is  undertaken 
to  obtain  the  solution  to  Equation  2-1,  This  approach  has  been  pursued  by 
other  investigators  to  solve  similar  problems  (see  References  2-4  to  2-7), 
Equation  2-1  can  be  witten  as  follows,  in  terms  of  four  linear  differential 
operators  of  the  first  order: 


°1°2°3V 


(2-2) 


where  (i*l,2,3,4)  denotes  the  linear  differential  operator 


,,3  3 

°i  “  3y  "  ^i  ax 


(2-3) 
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(2-4) 


and  are  the  roots  of  the  following  characteristic  equation: 

^ll^"^  -  -  2^26  ■■■  ^22  “  ° 

For  physically  meaningful  values  of  the  constants  solution  types  are 

possible: 

(1)  The  roots  of  Equation  2-4  are  two  pairs  of  complex  conjugates: 


yi  »  S-lg,  y2  «  Y+15,  y^  *  and  y^  *  y^ 
where  a  bar  denotes  a  complex  conjugate,  and  6, 6>0. 


(2-5) 


(2)  The  roots  are  pairwise  equal: 

hi  “  ^2  •  ^ 


(2-6) 


For  an  isotropic  plate,  y^^  ■  y2  •  1 


Rewrite  Equation  2-2  as  follows: 

D^F  -  gj-,  D3D^F  -  g^;  D2D3D^F  -  g^;  ^1  *  17  '  IT  “  ^ 
The  solution  to  the  last  equation  may  be  written  as: 
gi  >  F^(x+yiy) 


(2-7) 


(2-8) 


The  general  solution  for  F  may  then  be  written  as: 


(1)  In  the  case  of  different  complex  roots, 

F  -  Fj,  (x+y^y)  +  F2(x+y2y)  +  F3(x+yj^y)  +  F^(x+y2y) 


(2-9) 


(2)  In  the  case  of  pairwise  equal  complex  roots. 

F-  F, (x+y,y)  +  (x+y,y)  F,(x+y,y)  +>  (x+y,y)  +  (x+y,y) 

^  (2-10) 

Fi,  F2,  F3  and  are  arbitrary  functions  of  the  corresponding  variables. 
Redefining  these  variables  as: 

*1  -  *+Miy,  *2  “  *1  “  ^2  “  *■'■^2^  .  (2-11) 

and, recognizing  that  the  stress  function  I*  a  real  function  of  x  and  y. 
Equations  2-9  and  2-10  may  be  written  as; 


(1)  F  *  2Re  [fj(Zj^)  +  F2(z2)]  when  the  roots  are  different,  and 

(2)  F  =  2Re  [Fj^(Zj^)  +  Zj^F2(Zj^)]  when  the  roots  are  equal. 


(2-12) 


Any  expression  for  F  in  terms  of  arbitrarily  assumed  (z^^)  and 
$2  (z^)  functions  is  a  solution  to  the  cwo-dimensional  problem,  provided  these 
expressions  satisfy  the  appropriate  boundary  conditions.  Recall  that 
(;>j^  (zp  and  $2  governing  equation.  The  problem 

of  interest  involves  a  finite  anisotropic  plate  with  a  loaded  or  unloaded 
circuiar  or  elliptical  hole  (see  Figure  2-1).  In  general,  for  arbitrary 
geometry  and  boundary  conditions,  one  cannot  determine  closed  form  solutions 
for  (Jj  and  $2-  One  method  of  obtaining  approximate  solutions  is  to  consider 
series  expansions  of  the  functions  with  unknown  coefficients.  These  unknown 
coefficients  are  then  determined  by  satisfying  the  boundary  conditions  approxi¬ 
mately.  In  contrast  to  series  expansions  in  and  Z2,  faster  convergence  is 
obtained  if  series  expansions  are  assumed  in  coordinates  and  ^2*  obtained 
by  using  the  following  mapping  functions: 

-  (z^jt  yzj-a2-yjb2)/(a-iy^b)  »  1 

^2  “  ^  V  1  (2-18) 

In  a  laminate  with  a  circular  or  an  elliptical  hole,  the  Internal  boundary  gets 
transformed  to  an  ellipse  in  the  2^  ^2  (see  Equation  2-11),  The 

functions  in  Equation  2-18  map  the  Internal  boundary  to  a  unit  circle  in  the 
plane,  and  the  physical  region  of  the  laminate  to  the  exterior  of  the 
unit  circle  (see  Figure  2-1).  The  signs  of  the  square  root  terms  in  Equation 
2-18  are  chosen  such  that  the  Internal  boundary  is  mapped  on  to  the  unit  circle. 
Note  that  these  mapping  functions  are  analytic  functions,  and  hence  the  mapping 
is  conformal. 

In  the  5 j  ^2  P^®"®»  *^*2  assumed  to  be  the  following  modified 

Laurent  series  expansions: 


formal  Mapping  of  the  Finite  Anisotropic  Plate 
the  Zj-Z2  and  Cj  “  ^2 


30?St"5pp: 


V=i 


In  the  above  expressions,  and  are  coriplex  coefficients  v-hich  are  determined 
so  as  to  satisfy  the  boundary  conditions.  The  logarithmic  terms  drop  out  if  the 
force  resultants  on  the  internal  boundary  are  zero  (see  Reference  2-7). 


The  stresses 


ins  of  and  c^.and  these  derivatives  are  expressed 


as : 


4>,  -  (a  +  I  (na  C?  -  rta 

'  °  ft«l  /*  ' 


r-n»v,  J2  2  2^ 

))/  vz^-a  -y^b 


*2'  -"S.„?2"))/v4pI^ 


(2-20) 


In  specific  cases,  the  above  expressions  simplify.  For  example,  in  the  solu¬ 
tions  for  infinite  plates,  the  positive  exponent  terms  drop  out  as  the  stresses 
have  to  be  bounded  at  infinity.  For  most  infinite  plate  problems,  the  solu¬ 
tions  reduce  to  a  one  term  solution  (see  Reference  2-3). 

In  the  general  case  of  a  finite  anisotropic  plate  with  a  loaded  or  un¬ 
loaded  hole,  the  coefficients  can  be  numerically  calculated  to  satisfy  the 
boundary  conditions.  Generally,  a  convergent  solution  Is  obtained  by  truncating 
the  series  in  Equation  2-19  as  follows: 


In  order  to  determine  Che  finite  number  of  unknown  coefficients,  we  can 
select  the  same  total  number  of  points  on  the  inner  and  outer  boundaries,  and 
compute  the  unknown  coefficients  by  satisfying  the  boundary  conditions  exactly 
at  these  points.  However,  in  this  case,  the  calculated  solutions  at  other 
boundary  points  are  significantly  different  from  the  imposed  boundary  condi¬ 
tions  (see  Reference  2-4).  As  shown  in  Reference  2-4,  a  more  desirable 
approach  is  to  choose  a  larger  number  of  boundary  points  than  the  number  of 
unknown  coefficients,  and  to  satisfy  the  boundary  conditions  at  these  points 

12  "  '  ,  , 


in  a  least  squares  sense.  This  solution  procedure  has  been  adopted  in  the 
analysis  In  Reference  2-1.  Prior  tc  applying  the  least  squares  boundary 
collation  procedure,  the  single-valuedness  ot  the  displacements  and  the  rigid 
body  rotational  constraint  have  to  be  imposed. 

2. 3  Single-Valuedness  of  Displacem.ents 

In  the  series  expansions  for  and  (Equation  2-21),  the  logarithmic 
terms  are  multi-valued  and  give  rise  to  nulti-valued  displacements.  To  be 
physically  meaningful,  the  displacement  functions  have  to  be  single-valued, 
and  tnese  conditions  have  to  be  explicitly  imposed  on  the  coefficients  of  the 

logarithmic  terms. 

From  Equations  2-16  and  2-21,  the  contributions  of  the  logarithmic  terms 
to  the  displacements  u  and  v  are  seen  to  be  2Re(pja^ln^j^  +  and 

2Re (q^a^ln^j^  +  q26^1nC2).  respectively.  For  the  displacements  to  be  single¬ 
valued,  the  Increment  acquired  by  the  above  functions  in  describing  the  internal 
boundary  should  be  zero  (see  Reference  2-8),  On  the  internal  boundary, 

I  =  U2I  =1.  and  Cj  ^2  Therefore,  Re  |p^a^(ie)  +  P26jj(i0)  ^  -0 

and  Re  I  q^a^(i9)  +  q26Q(l9)  |  *=0,  Or, 

!m(q^ag  +  q^B^)  ■  q^a^  ♦  q^B^  -  q^a^ 

Imposition  of  the  above  two  conditions  will  ensure  that  the  computed  displace¬ 
ments  are  single-valued.  Equations  2-22  eliminate  two  of  the  unknown  constants: 


-  “ 

-  qjBj  ■  0  (2-22) 


“o'  l<'’2’rV)>“o  * 

«0  ■  *  <<llP2-’2'’l>\]/('llf>2-<l2h> 
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Zero  Rigid  Body  Rotation 


The  rotation  oj  of  an  infinitesimal  element  at  any  point  (x,  y)  in  the 
laminate  is  given  by  : 

“  3x  "  3y  (2-25; 

Using  Equation  2-16, 

w(x,y)  -  Re  ^U2P2  “ 

-  Re  |(UjPj-qpa^  [z^/  '4j-a^-Pjb^  +  l/Ca-iy^b)] 

+  (M2P2-‘l2^®l  ^^2^  } 

+  terras  involving  and  ^2  with  negative  and  positive  exponents  ^  2 

-  Re  [(yj^pj-qpaj/(a-lVjb)  +  (U2P2  -  q2)&i/(a-iU2b)]  (2-26) 

+  terras  which  are  functions  of  x  and  y  (i.e.,  are  not  constant) 

For  the  rigid  body  rotation  of  the  laminate  to  be  zero,  the  constant  terra 
above  should  be  zero.  Hence, 

Re  [{u^P^-q^)o,/{a-iwjb)  ♦  (y^Pg  -  q2)e/U-1P2‘>)]  -0 

Or,  (p^P^-q^)a^/{«-lM^b)  ♦ 

♦  (M^P^^q^)oj/(a  +  ♦  (W2P2*R2^^1^^*‘*^^^2^^ 

This  eliminates  one  more  unknown  constant: 


(2-27) 


2.5 


Least  Squares  Colutlon  to  an  Over-Determined  System  of  Linear  Equations 

Consider  the  following  ovei -determined  system  of  linear  algebraic 
equations: 

[a]  {X}  =  (b)  (2-29) 


where  [a]  Is  M  x  N,  (x)  is  N  x  1,  and  (b)  is  M  x  1  In  size,  and  M  >  N. 


The  least  squares  solution,  l.e.,  the  solution  that  minimizes  the  sum 


of  the  squares  of  the  error  in  each  equation,  is  given  by: 


I  ((  r  A..x.).bj")  « 
i-l  j-1  J  1  / 


0  for 


(2-30) 


Or, 


M  N 

r  (  I  A.  .x.-b.)A..  ■  0  fork -1,2 . N  (2-31) 

j«i  J  1  Ik 


The  above  N  equations  can  be  combined  to  produce  the  following  matrix  equation: 
[A]  {x}  -  [a]^  {b}  «  {0}  (2-32) 

The  least  squares  solution  is  then  expressed  as: 

{x}  -  |[a]^  ^A]|"^  [a]^  (b)  (2-33) 

2. 6  Least  Squares  Boundary  Collocation  Solution  Procedure 

Consider  the  finite  anisotropic  plate  in  Figure  2-1.  If  it  has  an  un¬ 
loaded  hole,  every  point  on  the  hole  boundary  is  stress-free.  That  is,  the 

normal  and  tangential  stress  components  (a_  and  T  «)  are  zero  along  the 

r  ry 

boundary.  In  this  case,  the  applied  loads  on  the  external  boundaries  are 
self-equilibrating.  Let  a  total  of  M  collocation  points  be  selected  on  the 
hole  boundary  and  the  external  boundaries.  At  each  point,  the  normal  and 
shear  stress  components  are  known.  Using  Equation  2-16,  all  the  2M  boundary 
conditions  are  written  in  terms  of  the  unknown  constants  in  Equation  2-21. 
Stress  boundary  conditions  on  the  hole  boundary  are  transformed  from  the  r-0 
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coordinates  in  which  they  are  easily  specified,  to  the  x  -  y  coordinates  in 
which  Equation  2-16  is  expressed,  using  appropriate  transformation  equations 
(see  Reference  2-2).  The  2M  boundary  conditions  are  expressed  in  matrix  form 


2Re.|  [B]  >  {($})  = 
or.  ;  [B  B]  J 


(2-34) 

(2-35) 


where  [b]  is  a  2M  x  (4N+2)  complex  matrix,  {6}  is  a  (4N+2)  x  1  vector  of 

the  unknown  complex  coefficients  (tt_jj^  . . .  ,6^,  •  .6j^)  , 

and  {R}  is  a  2M  x  1  real  vector  that  contains  the  specified  boundary  values. 

displacement  single-valuedness  and  tero  rigid  body 
rotation  conditions  (Equations  2-23,  2-24  and  2-28)  eliminates  three  unknowns. 
Using  the  foirm  in  Equation  2-35,  the  vector  of  unknown  complex  coefficients 
reduces  in  size  from  (8N+4)  x  1  to  (8N+1)  x  1.  The  reduced  system  of  equations 
is  written  as; 

tC]  W  =  (R)  (2-36) 

where  [c]  is  a :2M  x  (8N+1)  complex  matrix  and  {7}  is  the  vector  of  8N+1  un¬ 
known  complex  coefficients  °N*  ®  N’  ^-N’ . Oq. 

Oj  and  8^)  .  Rre-multiplying  equation  2-36  by  [c]^,  its  least  squares  solution  is: 


(7)  -  ([q^lcl)-'  [C]^{R) 


(2-37) 


Once  the  complex  coefficients  in  Equation  2-21  are  determined,  the 
stresses  and  the  displacements  are  calculated  using  Equation  2-16.  The 
accuracy  of  the  solution  is  determined  by  recalculating  the  stresses  at  the 
boundaries  and  comparing  them  to  the  imposed  boundary  conditions.  As  dis¬ 
cussed  la  ter  «  an  N  value  of  7  and  approximately  100  points  on  the  boundary 
are  sufficient  to  recover  the  imposed  boundary  conditions  with  5%.  , 


If  the  ahisotroplc  plate  in  Figure  2-1  has  a  loaded  hole,  the  assumptions 
made  in  Reference  2-9  ate  retained.  The  fastener  is  assumed  to  be  frictionless 
and  is  assumed  to  bear  over  half  the  hole  boundary.  The  fastener/laminate 
contact  problem  is  by-passed,  and  the  contact  solution  is  assumed  to  be  a 
cosinusoidal  distribution  of  the  radial  stress  (o^)  around  the  hole  (see 
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Figure  1-2).  The  tangential  stress  (^re>  is  zero  around  the  frictionless  hole 
boundary.  Results  from  recent  investigations  (References  2-5  and  2-6)  indicate 
that  the  contact  problem  could  affect  the  local  stresses  significantly.  Never¬ 
theless,  once  the  contact  stress  conditions  around  the  hole  boundary  are  com¬ 
puted  and  incorporated  into  the  least  squares  boundary  collocation  solution 
procedure,  the  two-dimensional  solution  is  computed  as  explained  earlier. 

The  discussed  boundary  collocation  solution  procedure  has  been  pro¬ 
grammed  to  be  the  FIGEOM  (Finite  Geometry)  computer  code  (see  Reference  2-1). 

2. 7  Effect  of  Number  of  Terms  in  the  Assumed  Series 

The  FIGECM  computer  code  was  Initially  used  to  determine  the  effect  of 

the  number  of  terms  in  the  assumed  series  expressions  (N  in  Equation  2-21)  on 

the  computed  solution.  For  this  purpose,  a  50/40/10  laminate  (from  Reference 

2-10)  with  a  1/4  inch  diameter  unloaded  hole  was  considered.  The  laminate  was 

subjected  to  a  tensile  loading  as  shown  in  Figure  2-2,  and  a  (0,D/2)/a  was 

X  o 

computed  for  various  values  of  N.  The  width  of  the  laminate  was  assumed  to  be 
small  (W/D=2)  to  Influence  the  computed  (0,D/2)  value,  A  plot  of 
(0,D/2)/a^  versus  N  indicates  that  tne  solution  coverages  when  N  exceeds  6  in 
value  (see  Fieure  2-2). 

Similar  studies  were  also  conducted  on  other  laminate  layups  and  on 
metallic  plates.  For  large  W/D  and  E/D  values  (outer  boundaries  located  with 
respect  to  the  center  of  the  hole  by  distances  in  excess  of  4D) ,  a  smaller  N 
value  (N  <  6)  yields  a  converged  solution.  Based  on  these  studies,  N  was  set 
equal  to  7  for  subsequent  analyses,  to  ensure  converged  solutions  under  all 
situations. 

2.8  Effect  of  Number  of  Collocation  Points 

Let  Njj  be  the  number  of  collocation  points  along  the  hole  boundary, 

and  let  Nj^  be  the  number  of  collocation  points  on  the  remote  boundaries.  The 

total  number  of  collocation  points,  M,  is  then  equal  to  N„  +  N^.  Referring 

n  R 

to  Section  2.6,  M  has  to  be  greater  than  4N  +  2,  where  N  is  the  number  of  terms 
in  the  assumed  stress  functions  (see  Equation  2-21).  When  N  ■  7,  M  has  to  be 
greater  than  30. 


50/40/10  Laminate  (h -0.1 21n.) 
W/D  -  2 

Open,  1/4  In.  die.  hole 


Number  of  Terms  in  the  ^  Series,  N 


Figure  2~2.  Effect  of  N  on  Solution  Convergence 


Con?  Icr  a  50/A0/10  laminate  subjected  to  a  static  tensile  load  in 

double  sheax  (see  Figure  2-3).  The  applied  load  is  fully  resisted  by  the 

fastener.  Three  situations  are  analyzed  to  evaluate  the  effect  of  the  number 

of  collocation  points  on  the  convergence  of  the  solution.  The  first  situation 

involves  a  5/16  inch  diameter  hole  and  a  laminate  geometry  defined  by  E/D  =  1.  2 

and  W/D  1.6.  The  second  situation  Involves  D  =  1/4  inch>  E/D  =  1.5  and  W/D  =  2.0. 

The  last  example  considers  a  laminate  with  D  =  1/4  inch,  E/D  =  3  and  W/D  =  4. 

In  every  case,  the  computed  value  at  (0,  D/2)  is  normalized  with  respect  to 

the  applied  remote  stress  value  (1  ksi)  to  obtain  a  stress  concentration  value 

of  K.  The  results  for  the  three  laminate  geometries  are  presented  in  Figure  2-3. 

Every  remote  boundary  contains  N  /4  collocation  points.  It  is  seen  that,  when 

E/D  and  W/D  are  very  small  (1.2  and  1.6,  respectively),  K  converges  slowly  with 

N  .  But,  when  E/D  and  W/D  are  not  very  small  (E/D  i  1.5,  W/D ^2),  K  converges 
K 

more  rapidly  when  is  Increased.  N„  =  30  and  N„  “  50  provide  approximately 
K  H  n 

the  same  results  (first  example).  If  E  ^  1.5  D  and  W  ^  2D,  K  changes  by  less 

than  5%  when  N„/4  is  Increased  beyond  a  value  of  10. 

K 

rased  on  the  results  in  Figure  2-3,  subsequent  analyses  using  the 

FIGEOI  computer  code  were  carried  out  with  N  •  50  and  N_/4  »  10,  using  a 

H  R 

total  of  90  collocation  points,  unless  specified  otherwise. 

2.9  Sample  FIGEOM  Predictions 

Having  established  N,  N  and  N  values  (7,  50  and  40,  respectively) 

n  R 

for  solution  convergence,  FIGEOM  Is  now  applied  to  a  few  sample  test  cases 
to  demonstrate  its  predictive  capability.  The  first  example  considers  a 
50/40/10  laminate  with  an  unloaded  hole  (Reference  2-10),  subjected  to  a 
uniform  tensile  loading.  The  hole  diameter  Is  1/4  inch  and  E/D  and  W/D  are 
small  (2).  Figure  2-4  compares  the  normalized  (0,  y)  variation  across 
the  hole,  predicted  by  FIGEOM,  with  the  predictions  using  two  other  analyses 
(References  2-7  and  2-9).  FIGEOM  predictions  agree  well  with  those  based  on 
a  similar  analysis  (Reference  2-7),  but  a  considerable  difference  Is  observed 
In  comparison  to  the  infinite  plate  solution  (Reference  2-9).  The  infinite 
plate  solution  will  be  approximately  equal  to  the  finite  plate  solution 
(FIGEOM)  when  E/D  and  W/D  values  are  much  larger  than  four. 
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(0,D/2)  /a 


50/40/10  Laminate  (Ref.  2-10);  h  «  0.119  in. 

L  ■Sin.;  Full  bearing  in  double  shear  (0,15  in.  Al) 

N  ••  number  of  collocation  points  on  the  remote  boundaries 
Nu  “  number  of  collocation  points  on  the  hole  boundary 


Effect  of  Number  of  Collocation  Points  on 
Solution  Convergence. 


The  effect  of  plate  width,  or  the  W/D  ratio,  on  the  stress  concentration 
at  (0,  D/2)  is  shown  in  Figure  2-5.  The  unloaded  hole  is  1/4  inch  in  diameter 
and  E/D  is  equal  to  3.  Three  laminates  from  Reference  2-10  are  considered. 

It  is  seen  that  the  stress  concentration  is  relatively  unchanged  beyond  W/D= 4, 
but  increases  significantly  when  W/D  is  decreased  below  four. 


The  effect  of  the  edge  distance  on  the  stress  concentration  l.':.  (0,  D/2) 
is  shown  in  Figure  2-6.  In  this  case,  E/D  has  relatively  no  effect  on  the 
stress  concentration  in  all  three  laminates.  If  the  hole  had  been  a  loaded 
hole,  E/D  would  have  had  a  more  significant  effect  on  K. 

The  effects  of  E/D  and  W/d  on  the  tangential  stress  distribution 
around  the  boundary  of  a  loaded  hole  are  shown  in  Figures  2-7  and  2-8.  The 
solutions  from  Reference  2-9,  in  both  the  figures,  correspond  to  a  location 
0.02  inch  from  the  hole  boundary,  and  therefore  yield  lower  stress  concentra¬ 
tions.  The  finite  plate  solution  (FIGEai)  yields  a  larger  stress  concentra¬ 
tion  at  6  ■  90°  in  compariscn  to  the  infinite  plate  solution,  the  approximate 
finite  width  solution  and  the  finite  element  solution  presented  in  Reference 
2-9.  FIGEOM  predictions,  corresponding  to  a  location  0.02  inch  away  from  the 
hole  boundary,  are  expected  to  agree  well  with  the  finite  element  solutions. 

At  9  »  90°,  Og  ■  (Oi  D/2).  If  E/D  is  reduced  below  three  and  W/D  is  re¬ 
duced  below  six,  the  effect  of  the  closer  boundaries  will  result  in  much 
larger  stress  concentrations  at  0  =  90°. 

^  plate  geometry  is  shown  in  Figure  2-9.  This  sample 

considers  a  laminated  lug  that  transfers  the  applied  load  to  aluminum  plates 
in  a  double  shear  configuration.  The  lug  has  a  50/40/10  layup,  and  the  effect 
of  R/ r  on  its  stress  concentration  at  x»  o  and  y  *  r  is  shown  in  Figure  2-9. 

In  this  case, N/4  was  selected  to  be  20  Instead  of  the  value  (10)  used  in  other 
Sample  analyses.  The  results  in  Figure  2-9  indicate  that  a^(o,r),  normalized 
with  respect  to  the  average  bearing  stress,  decreases  when  the  outer  boundary 
is  moved  away  from  the  loaded  hole  (R/ r  increases) .  A  rectangular  geometry 
at  R/r>>  2  yields  a  lower  normalized  stress  than  the  non-rectangular  geometry. 

The  final  example  considers  a  situation  where  a  fraction  of  the  total 
applied  tensile  load  is  transferred  directly  to  the  fastener  hole  location. 
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70/20/10 


bearino 


bearing 


0  -  0.25  in. 

E/D  >  9 
W/0  >  8 

60/40/10  Layup 


At  the  hole  boundary 


At  a  location  0.02 
in;  from  the  hole 
boundary  — — ^ ^ 


Finite  Element  Solution  (Ref.  2-8) 
infinite  Plate  Solution  (Ref.  2-9) 
FIQEOM 


?20  40  60  80  100  120  140  160  180 

CIRCUMFERENTIAL  LOCATICN  AROUND  THE  HOLE  BOUNDARY. 

9  (DEGREES) 


Figure  2-8.  Effect  of  E/D  ■  9  and  W/D 
with  0. 


8  on  the  Oq  Variation 


(0,  D/2)  /o 


50/40/10  laminate;  h  =  0.119  in.;  double  shear  load 

transfer  to  0. 155  in.  A1  plates;  P/Wli  =  1  ksi 

D  “  1/4  in.;  a,  =  P/Dh  =  W/D  ksi 
brg 


Figure  2-9.  Effect  of  Lug  Geometry  on  Stress  Concentration, 
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The  difference  between  the  total  applied  load  and  the  fastener  load  is  referred 
to  as  the  by-pass  load.  This  example  represents  an  isolated  fastener  location 
in  a  plate  that  is  mechanically  fastened  to  another  using  many  bolts  (see 
Figure  1-1).  Results  for  three  laminates,  with  D  =1/4  inch,  E/D  =  5  and 
W/d  •  8,  are  presented  in  Figure  2-10.  Presented  results  indicate  that  ':he 
stress  concentration  increases  when  the  fractional  bolt  bearing  load  is 
increased.  The  predicted  linear  increase  in  K  with  an  increase  in  the  ir actional 
bolt  load  provides  an  explanation  for  the  linear  strength  reduction  with  the  bolt 
load  fraction  observed  in  Reference  2-10. 

2.10  Summary 

A  two-dimensional  analysis  (FIGEtti)  for  finite  anisotropic  plates  with 
loaded  or  unloaded  holes,  developed  in  Reference  2-1,  was  discussed.  The 
analysis  uses  a  boundary  collocation  technique  and  computes  the  solution  in  a 
least  squares  sense.  Examples  were  presented  to  demonstrate  the  capability  of 
the  analysis  in  computing  the  effect  of  plate  geometry  on  the  stress  concentra¬ 
tion  at  the  boundary  of  the  loaded  or  unloaded  holes.  Computed  solutions  were 
compared  with  Infinite  plate  solutions  to,  demonstrate  the  significant  Increase 
in  the  stress  .concentration  when  the  outer  plate  boundaries  are  moved  closer 
to  the  hole.  ' 

The  major  limitation  of  the  FIGEOM  computer  code  is  its  approximation 
of  the  fastener/plate  contact  stress  distribution.  Though  the  contact  stress 
distribution  and  the  contact  region  depend  on  the  laminate  layup,  FIGEOM 
assumes  a  cosinusoidal  bearing  stress  distribution  over  half  the  hole  boundary. 
Nevertheless,  the  developed  analysis  is  a  significant  improvement  over  the 
infinite  plate  analysis  presented  in  Reference  2-9. 
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SECTION  3 


FASTENER  ANALYSIS 

3.1  Introduction 

In  the  computation  of  the  two  dimensional  stress  state  in  a  bolted  plate, 
it  Is  generally  assumed  that  the  fastener  Is  rigid  and  that  the  fastener 
plate  displacement  due  to  Inplane  loads  does  not  vary  In  the  plate  thickness 
direction  (see  Section  2) .  In  most  of  the  practical  situations,  this 
assumption  is  not  valid,  and  the  stress  field  at  the  fastener  location  Is 
complex  and  three  dimensional  In  nature.  The  three  dimensional  stress  field 
at  the  fastener  location  ^  a  bolted  metal  or  composite  plate  Is  influenced  by 
many  factors:  (1)  fastener  size,  (2)  fastener  stiffness,  (3)  fastener  end  con¬ 
straints;  (4)  fastener  torque,  (S)  hole  clearance,  (6)  properties  of  the 
bolted  plates,  (7)  stacking  sequence  of  the  bolted  laminate,  (8)  load 
eccentricity  Induced  by  Joint  configuration,  etc. 

V  3-1  and  3-2  Illustrate  the  difference  between  a  double  lap  and 

a  single  lap  joint  configuration.  The  single  lap  configuration,  due  to  the 
eccentricity  4n  the  load  path,  will  affect  the  local  stress  field  more 
significantly.  This  effect  is  further  influenced  by  the  fastener  properties. 

If  the  fastener  modulus  ^  large  compared  to  the  bolted  plate  modulus,  and 
the  fastener  diameter  is  large  compared  to  the  plate  thickness,  the  fastener 
bending  and  shear  stiffnesses  will  be  large  enough  to  cause  it  to  act  like  a 
rigid  fastener.  Otherwise,  fastener  bending  and  shear  deformation  will 
influence  the  local  stress  field  significantly  (see  Figure  3-3). 

Fastener  end  constraints  also  have  a  significant  effect  on  the  local  stress 
state.  Figure  3-4  shows  three  situations  of  interest.  The  protridjng  head 
fastener  with  a  large  applied  torque  value  (Figure  3-4a)  essentially  creates 
a  nearly  fixed  end  boundary  condition  (constraint).  A  pin  permits  free 
rotation  at  the  boundary  (Figure  3-4c),  and  a  highly  torqued  countersunk 
fastener  creates  nearly  fixed  and  nearly  free  constraints  at  the  nut  and 
head  locations,  respectively  (Figure  3-4b).  Intermediate  torque  values 
can  be  represented  by  elastic  constraint  equations  that  quantify  constraints 
between  fixed  and  free  conditions. 
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A.  V 


Figure  3-2.  A  Single  Lap  Configuration 
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El  and  AG  are  comparable 
to  plate  EA 


(b)  Flexible  Fastener —  Measurable  fastener  bending 

and  shear  deformation 


EA  -  plate  axial  stiffness 
EX  -  fastener  bending  stiffness 


GA  -  fastener  shear  stiffness 


Figure  3-3 •  Typical  Rigid  and  Flexible  Fasteners. 


El,  AG  are  large  compared  to 
plate  EA  value 


D/h  »1 


(a)  Rigid  Fastener  —  negligible  fastener  bending  & 

shear  deformation 


(a)  Fixed  -  Fixed  Conditions  (Protruding  head  fastener  — 

high  torque-up) 


(b)  Fixed  -  Free  conditions  (countersunk  fastener,  —  high 

torque-up) 


(c)  Free-Free  Conditions  (pin) 

Figure  3-4.  A  Single  Lap  Configuration  with  Various 
End  Constraints  on  the  Fastener. 
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Hole  clearance  effects  are  also  obvious  from  Figure  3-4.  In  the  case  of 
a  free  pin,  for  example,  an  Initial  hole  clearance  will  cause  the  pin  to 
rotate  in  place  or  "cock  "  when  a  lord  is  applied.  If  a  protruding  head 
fastener  is  used,  the  friction  due  tc  :he  applied  torque  is  overcome  prior 
to  fastener  cocking  (rigid  body  rotate  an). 

fastener  bending  and  shea  '  deformation  effects  are  not  negligible 
or  if  fastener  rotation  is  Induced  by  any  of  the  discussed  factors,  the  local 
Btress  field  will  be  affected  by  the  stacking  sequence  of  the  bolted  laminate. 

This  is  due  to  the  non-uniform  strain  distribution  in  the  thickness  direction 
of  the  bolted  plate. 

^  following  subsections  present  details  of  an  analysis  that  was  developed 

to  predict  the  three  dimensional  stress  state  at  a  fastener  location  in  a 
bolted  laminated  or  metallic  plate.  The- discussed  analysis  is  similar  in 
approach  to  that  in  References  3-1  to  3-3.  But,  while  the  analysis  in 
References  3-1  to  3-3  is  restricted  to  bolted  metallic  plates,  the  analysis 
developed  in  this  program  is  applicable  to  metallic  and  laminated  anisotropic 
plates. 

3.2  Summary  of  Approach 

A  brief  summary  of  the  analytical  approach  is  presented  with  the  aid  of 
Figure  3-5.  For  a  single  lap  joint  configuration  subjected  to  a  tensile  load 
(Figure  3- 5a)y a  typical  fastener/bolted  plate  displacement  variation  is 
shown  in  Figure  3-5(b) i  The  distribution  of  the  contact  force  is  influenced 
by  the  many  factors  discussed  in  Section  3.1,  and  is  not  a  continuous 
function  along  the  fastener  axial  direction  (or  the  plrte  thickness  direction) 
for  a  laminated  bolted  plate.  Figure  3-5(b)  shows  a  typical  contact  force 
distribution  in  bolted  metallic  plates.  The  resultant  of  the  contact  force 
distribution  will  be  equal  to  the  applied  load  (P)  in  magnitude,  but  will 
not,  in  general,  lie  along  the  line  of  action  of  P.  This  is  because  its  line 
of  action  is  determined  based  on  moment  equilibrium  considerations.  A  free 
body  diagram  of  the  fastener  will  include  contact  forces  that  are  opposite 
in  sense  to  those  shown  in  Figure  3-5(c).  The  spring  constants  represent  the 
resistance  offered  by  the  bolted  plate  to  fastener  displacement.  In  a 
laminated  plate,  spring  constmts  vary  from  ply  to  ply  and  are  dependent  on  ply 
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(a)  Single  Lap  Bolted  Joint 


(b)  Typical  Fastener/Plate  Displacement 
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Figure  3-5. 


Mathematical  Representation 

Representation  of  a  Single  Lap  Configuration 
by  an  Equivalent  Fastener  Problem. 
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fibM  orientations.  In  a  metallic  plate,  the  spring  constant  will  be  invariant 
in  the  thickness  direction.  The  various  springs,  with  appropriate  constants 
assigned  to  each,  mathematlcaliy  replace  the  bolted  plates  by  an  elastic 
foundation  whose  modulus  id  piecewise  uniform  in  general. 

The  three  dimensional  stress  field  at  the  fastener  location  is,  therefore, 
computed  by  obtaining  the  solution  to  a  mathematical  problem  that  represents 
the  fastener  as  a  beam  resting  on  an  elastic  foundation  with  piecewise 
uniform  moduli  that  tepresait  the  various  plies  in  a  bolted  laminate.  The 
fastener  is  modeled  as  a  Timoshenko  beam  to  account  for  bending  and  shear 
deformation  effects.  The  foundation  is  represented  by  a  general  bilinear 
contact  load  versus  deflection  curve,  to  account  for  a  reduced  ply  stiffness 
after  an  initial  damage  (local  failure)  is  precipitated.  At  the  head  and 
nut  locations  of  the  fastener,  rotational  constraints  are  Introduced.  These 
constraints  are  influenced  by  the  applied  torque  and  the  size  of  the  washers, 
if  any.  Fastener  torque  also  introduces  friction  forces  on  the  bolted  plates, 
reducing  the  load  transferred  directly  by  the  fastener.  If  either  of  the 
bolted  plates  is  a  laminate,  its  stacking  sequence  will  Influence  the 
fastener  displacement,  and  hence  the  load  distribution  in  the  thickness 
direction. 

The  governing  differential  equation  for  the  fastener  displacement  is 
solved  after  the  derivatives  are  replaced  by  appropriate  differences.  The 
finite  difference  formulation  of  the  problem  facilitates  fast  and  accurate 
computations  via  an  easily  automated  solution  procedure.  For  a  symmetric 
double  lap  configuration,  symmetry  conditions  at  the  center  replace  the  nut 
location  boundary  conditions  used  otherwise.  Solutions  are  obtained  over 
the  region  of  contact  in  each  plate,  and  continuity  conditions  are  enforced 
where  the  two  plates  are  in  contacts  The  following  8ub->8ection8  present 
details  of  the  analysis. 

3.3  Governing  Differential  Equation 

The  fastener  is  modeled  as  a  Timoshenko  beam  to  account  for  shear 
deformation  effects.  Figure  3-' 6  shows  the  deformed  state  of  an  infinitesimal 
segment  of  the  fastener,  and  Figure  3-7  presents  the  assumed  sign  conventions 
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for  transverse  shear  (V)  and  bending  moment  (M)  in  the  fastener.  The  coordi¬ 
nates  along  the  fastener  axis  and  the  loading  direction  (in  the  plane  of  the 
bolted  plate)  are  labeled  z  and  x,  respectively.  Under  load,  a  plane  section 
AB  deforms  to  the  form  ^  shown  in  Figure  3-6.  In  doing  so,  the  cross- 
section  undergoes  a  translational  displacement  u(z)  in  the  x  direction,  a 
bending  rotation  ^(z),  and  shear  deformation.  While  t{)  is  not  a  function  of 
X,  the  shear  strain  due  to  y  is,  due  to  the  variation  of  the  transverse 
shear  stress  (t  )  in  the  x  direction;  This  results  in  the  curved  shape  for 

ZX 

AB.  Representing  the  average  cross-sectlonal  shear  distortion  (rotation)  by 
the  symbol  ()),  the- following  relationship  is  assumed  (Timoshenko  beam  theory): 

V  -  / 

where  A  is  the  fastener  cross-sectional  area,  G  is  its  shear  modulus,  and  X 
,  is  a  shear  correction  factor  accounting  for  nonlinear  distribution  in  the 
X  direction.  The  total  rotation  of  the  section  AB  is  denoted  by  u* ,  where  the 
prime  denotes  differentiation  with  respect  to  n.  From  Figure  3-6,  it  follows 
.‘-that: 

u*  ■  (j)  -  ij)  (3-2) 

The  bending  moment  at  any  z  location  is  expressed  as  follows: 

M  -  /  a  xdA  -  /  (E  e  )  xdA  -  /E(xifi*)  xdA  -  Elij;'  (3-3) 

z  z  .. 

where  a  is  the  fastener  bending  stress,  e  is  the  axial  strain  in  the  fastener 
due  to  bending,  E  is  the  fastener  Young's  modulus,  and  1  is  the  moment  of 
inertia  of  the  fastener  cross-section  about  the  y  axis  (normal  to  the  rz 
plane) . 

If  the  fastener  is  subjected  to  a  distributed  transverse  load  q(z),  force 
and  moment. equilibrium  considerations  yield  the  following  relationship  (see 
Figure  3-7) : 


■;  .38 


(3-4) 


V  =-q,  and 
m'  -  V 


(3-5) 


Again,  primes  denote  differentiation  with  respect  to  z.  Equations  3-1  to 
3-5  yield  the  following  relationships  for  M,  V  and  ili: 


M  -  -El  [u**  +  q/(XGA)] 

V  -  -El  [u"’  +  q'/(AGA)] 

ij;  -  -(EI/XGA)  [u"’-'qV(XG/)]  -u’ 


(3-6) 

(3-7) 

(3-8) 


Equations  3-4  and  3-5  may  be  combined  to  yield  the  ^  llowlr^  equilibrium 
equat ion: 


M  “  -q 


(3-9) 


Substituting  Equation  3-6  into  the  above  equation,  the  following  governing 
equation  la  obtained: 


u"”  +  q”  /(XGA)  -  q/(El)  -  0 


(3-10) 


This  equation,  where  q  ■  q(z),  governs  the  displacement  of  the  fastener, 
u""  ■  d^u/dz^  and  q"  -  d^q/dz^. 

3.4  Monlinear  Foundation  Behavior 

The  effect  of  the  bolted  plate  (metallic  or  laminated)  on  the  fastener 
displacement  is  represented  by  the  transverse  loading  term  q(z)  In  the 
Equation  3-10.  The  q(z)  terra  is,  In  turn,  linearly  related  to  the  fastener 
displacement  u(z)  through  the  foundation  modulus  k(z).  For  the  more  general 
laminated  foundation,  the  foundation  modulus  varies  from  ply  to  ply,  and  la 
uniform  within  a  ply.  k(z)  is,  therefore,  piecewise  uniform. 

In  the  developed  analysis,  every  ply  is  also  assumed  to  be  a  bilinear 
elastic  (Hencky)  material  (see  Figure  3-8).  This  representation  of  the  ply 
behavior  permits  the  prediction  of  a  local  damage  In  the  ply  (when  u  >  u^) 
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Figure  3-8.  Bilinear  Elastic  Behavior  of  a  Fly 


that  is  not  catastrophic.  The  ply  modulus  is  for  0  <  u  <  u^.  Depending 
on  the  type  of  damage  predicted  in  the  ply,  its  modulus  (Kt)  beyond  u  =  u  is 
set  to  be  greater  than,  equal  to  or  less  than  zero.  If  k2  ^  Ot  the  ply 
exhibits  a  hardening  behavior;  if  ^2  <  0,  it  exhibits  a  softening  behavior; 
and  if  k2  “  0,  the  ply  is  an  elastic-perf ectly  plastic  material.  Wnen  u 
takes  the  value  of  u*  (Figure  3-8),  the  ply  loses  its  load-carrying  capability. 
At  this  load  level,  the  ply, is  assumed  to  have  failed  totally  and  its  modulus 
and  load  are  reset  to  zero.  Adjacent  unfailed  plies  share  the  load  that  is 
released  by  the  totally  failed  ply. 

The  general  ply  load  versus  displacement  behavior  (Figure  3-8)  is 
expressed  mathematically  as; 


q(z)  *  -ku  -  ku 


(3-11) 


where,  k  ■  k^^  for  an  undamaged  ply 

k  ■  k^  ■  for  a  partically  damaged  ply 
k  ■  0  for  a  totally  damaged  ply 

k  -  kj^  -  k2  •  (l-a)k^  for  a  partically  damaged  ply 
k  ■  0  for  a  totally  damaged  ply 
u  «  0  for  an  undamaged  ply 

u  ■  for  a  partially  or  totally  damaged  ply 
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(k^,  u^,  k2  and  u*  or  Pyj^tj^te^^initlal^  fully  define  the  general  ply  behavior. 
The  computation  of  k^  for  the  various  plies  is  discussed  in  the  following 
sub-section,  u^  and  u*  are  dependent  on  the  failure  criteria  used  to 
predict  partial  and  total  ply  damage.  k2  *  0dCj|^  will  be  established  by 
assigning  a  values  for  the  various  partial  damage  types.  If  the  ply  behavior 
can  be  adequately  represented  by  a  linear  elastic  approximation,  a  simplified 
form  of  Equation  3-11  may  be  used. 


3.5  Computation  of  Initial  Foundation  Moduli  (kj^)  for  the  Various  Ply  Types 

If  a  bolted  plate  is  a  laminate,  the  various  fiber  orientations  in  its 
lay-up  determine  the  number  of  ply  types  in  tue  laminate,  assuming  all  plies 
tc  be  made  of  the  same  material.  A  metallic  bolted  plate  has  only  one  ply 
type.  The  Initial  foundation  modulus  (k^  in  Figure  3-8)  for  each  ply  type 
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Is  computed  by  considering  a  situation  where  the  local  three  dimensional 
stress  field  can  be  approximated  by  a  two  dimensional  stress  field.  This  is 
done  by  assuming  that  the  fastener  is  rigid  in  bending  and  in  shear 

(GA**^  “»),  and  that  the  load  is  transferred  in  a  symmetrical  double  lap  con¬ 
figuration.  In  this  situation,  the  inplane  strains  do  not  vary  in  the 
thickness  direction  when  only  inplane  loads  are  present.  This  strain  field 
and  the  corresponding  stress  state  in  each  ply  are  computed  using  the  FIGEOM 
computer  code'  (Ref erence  2-1) . 

An  effective  fastener  displacement  in  the  load  direction  (x)  is  obtained 
prior  to  computing  for  the  various  ply  types.  This  is  accomplished  by 
equating  the  total  work  done  by  the  fastener-imposed  radial  bearing  stresses 
around  the  hole  boundary  (ZW^)  to  the  expression  (1/2)  p  where  P  is 

the  total  fastener  load  and  is  the  effective  fastener  displacement  in 


a  hRAe  Cos(e.+  O.5A0  )]  (3-12) 
i«l  ■“  "  1-1  “  ^  ^  i 

where  n  is  the  number  of  collocation  points  around  the  hole  boundary,  and 
Figure  3-9  describes  the  various  quantities  that  are  present  in  Equation  3-12. 
R  is  the  radius  of  the  hole  and  h  Is  the  total  plate  thickness. 

The  strains  corresponding  to  a  fastener  load  of  P  are  invariant  in  the 
plate  thickness  direction.  Using  the  appropriate  stiffness  matrix  for  each 
ply  type  and  these  strains,  ply  stresses  are  computed  at  each  collocation 
point.  Stresses  in  rectangular  coordinates  (x,y)  are  transformed  to  stresses 
in  polar  coordinates  (R, 6)  using  appropriate  transformation  relationships.  The 
load  in  the  J  th  ply,  in  the  load  direction,  is  computed  as  follows 

(see  Figure  3-10); 


the  load  (x)  direction: 
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Figure  3-9.  Computation  of  Work  Done  by  the  Fastener  Bearing 
Stress 
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Figure  3-10.  Computation  of  Ply  Load. 
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(3-13) 


The  Initial  foundation  modulus  (kj^)  for  the  j  th  ply  type  is  then  computed 
using  the  followlngvreJ.'tionship; 


^xj'^^^j^eff^ 


(3-14) 


where  and  are  defined  in  Equations  3-12  and  3-13,  respectively,  and 
hj  is  the  thickness  of  the  j  th  ply  type. 

3.6  Boundary  and  Continuity  Conditions 

The  boundary  and  continuity  conditions  on  a  fastener  that  bolts  two 
plates  in  a  single  lap  configuration  are  shown  in  Figure  3-11.  The  portion 
of  the  fastener  in  eacit  plate  is  shown  separately.  The  load  (P)  in  each 
plate  is  enforced  as  a  shear  boundary  condition  at  the  interfacial  locction, 
to  satisfy  force  equilibrium  requirement.:  The  shear  force  values  at  ti  i 
outer  boundaries  (the  head  and  nut  locations  of  the  fastener)  are  set  equal 
to  zero,  since  the  load  transfer  is  effected  between  these  locations. 

At  the  interface  between  the  "bblted  plates,  continuity  of  the  bending  slope 
and  the  bending  moment  are  enforced.  Continuity  of  displacement  is  not 
enforced  at  this  location. >  This  is  because  u  (z)  represents  the  fastener/ 
plate  displacement,  and  undergoes  a  finite  discontinuity  across  this  inter- 
facial  location  in  the  Joint. 

At  the  fastener  head  location,  the  head  type  and  the  presence  of 
washers,  if  any.  Influence  the  constraint  against  free  rotation.  A  washer 
and  a  nut  offer  a  similar  constraint  at  the  other  boundary.  The  constraints 
at  the  fastener  head  and  nut  locations  can  be  generalized  as  shown  in 
Figure  3-11,  where  and  R2  quantify  the  elastic  restraint.  For  a  pin- 
connected  Joint  (Figure  3-4c),  for  example,  R^  and  R2  are  set  equal  to  zero. 


For  a  countersunk  fastener,  is  set  equal  to  zero  (at  the  head  location). 
Fastener  torque  influences  and  R2  values  significantly.  If  a  protruding 

head  fastener  is  hi^ly  torqued  (Figure  3-4a),  R^^  and  R2  values  become  very 
large  (’^  ®  ),  and  tend  to  reduce  the  bending  slope  to  zero  at  the  head  and 
nut  locations* 

^  elasticity  analysis  that  could  be  useful  in 

the  determination  of  the  rotational  constant  R  in  M  >  RiJ^.  The  analysis 
assumes  a  rigid  circular  disk  to  be  supported  by  a  transversely  Isotropic 
half  space,  and  subjected  to  a  transverse  load  and  a  bending  moment.  The 
derived  expression  for  R  is: 

R  =  8Gr^/(3«{)^)  (3-15) 

where  G  is  the  transverse  shear  modulus  of  the  half  space,  r^  is  the  radius 
of  the  disk,  and  4)^  is  a  dimensionless  constant  that  varies  from  0.7  for 
an  isotropic  space  to  1.0' for  a  highly  anisotropic  space.  In  a  practical 
' joint,  an  annular  disk  (washer)  is  supported  by  a  foundation  of  finite 
thickness  and  the  expression  in  equation  3-15  is  not  directly  applicable. 

An  empirical  expression  will,  therefore,  be  derived  for  R,  based  on  a 
correlation  between  analysis  and  experiment.  This  expression  will  relate 
R  to  the  washer  diameter  and  the  fastener  torque. 

Load  transfer  in  a  symmetric  double  lap  configuration  yields  the  boun¬ 
dary,  symmetry  and  continuity  conditions  shown  in  Figure  3-12.  The 
boundary  conditions  at  the  fastener  head  location,  and  the  shear  and  con¬ 
tinuity  conditions  at  the  interface  between  adjacent  plates,  are  identical 
to  those  discussed  earlier  (see  Figure  3-11).  The  assumption  of  a 
syssmetric  double  shear  situation  requires  the  imposition  of  the  symmetry 
conditions  shown  in  Figure  3-12.  Symmetry  Introduces  a  zero  shear 
a  zero  bending  slope  condition  at  the  midplane  of  plate  2.  It  is  noted 
that  only  half  the  fastener  is  analyzed  for  a  symmetric  double  shear 
situation. 

A  total  of  8  boundary/ continuity/ symmetry  conditions  are  identified 
for  each  joint  configuration  —  the  single  lap  joint  in  Figure  3-11,  and 
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Figure  3-12.  Boundary  and  Continuity  Conditions  for  a 
Typical  Double  Lap  Joint. 
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the  double  lap  joint  in  Figure  3-12.  For  the  single  lap  joint  (Figure  3-11), 
the  8  conditions  are; 

V  =■  0  and  M  =■  at  the  fastener  head  location 

V  ■  -P  at  the  plate  1/plate  2  Interface,  in  the  plate  1  region  V 

V  «  -F  at  the  plate  l/plate  2:  interface,  in  the  plate  2  region  I 

M  at  the  plate  l/plate  2  interface  in  the  plate  1  region  ■>  M  at  the/ 

plate  1/plate  2; Interface  in  the  plate  2  region  )(3-16) 

at  the  plate  l/plate  2  interface  in  the  plate  1  region  ■<>  at  thei 

plate  l/plate  2  Interface  in  the  plate  2  region  I 

V  =  0  aT'H  M  “  at  the  nut  location  / 

For  the  double  lap  joint  (Figure  3-12),  the  first  6  conditions  are  the 
same  as  above.  The  last  two  conditions  at  the  nut  location  are  replaced  by 
the  following  conditions  at  the  midplane  of  plate  2*. 

V  •>  0  and  M  >  0  (4'  ■  0)  at  the  midplane  of  plate  2  (plane  of  symmetry) 

(3-17) 

3.7  Finite  Difference  Formulation  of  the  Problem 

The  governing  equation  for  the  fastener/plate  displacement  (Equation 
,3-10)  containsfexpresslons  that  are  dependent  on  the  relationship  between 
'.the  ply  load  and;  the  ply  displaces  3-11).  Incorporation 

of  Equation  (3-11);  into  Equation  (3-10)  will  result  in  a  fourth  order, 
ordinary  differential  equation  for  u(z)  with  variable  coefficients 
(because  k  isi  a  function  of  z).  A  finite  difference  approximation  of  the 
governing  equation'  (and  the  boundary  and  continuity  conditions)  is  adopted 
to  obtain  the  fastener/plate  displacement  and  the  corresponding  fastener 
load  distribution  in  the  thldkness  direction  of  the  bolted  plates.  This 
provides  a  rapidly:  executable  solution  scheme  that  can  be  economically 
executed  many  times  to  predict  progressive  failures  in  bolted  joints. 

Details  of  the  finite  difference  formulation  of  the  problem  are  presented 
below. 


Figure  3-13.  An  Example  of  the  Node  Layout  and  Numbering 

Scheme  in  a  Single  Lap  Shear  Joint  Configuration. 
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Figure  3-lA.  A  General  Node  Arrangement  with  n  Nodes  in  Plate  1 
and  m  Nodes  in  Plate  2. 
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Figure  3-13  describes  how  the  continuous  fastener  is  discretized 
-o  a  finite  number  of  nodal  points.  In  the  chosen  example,  plates  1  and 
.  are  assumed  to  be  5-ply  and  4-ply  laminates,  respectively,  bolted  together 
in  a  single  lap  configuration.  The  portion  of  the  fastener  in  each  plate 
is  represented  by  equally-spaced  nodes  located  at  the  ply  midplanes.  If 
either  plate  is  metallic,  it  Is  divided  into  m  or  n  number  of  plies  of  the 
same  properties.  To  enable  the  use  of  a  central  difference  scheme  for 
thfe  governing  equation  and  the  boundary  conditions,  two  'false'  nodes  are 
assumed  to  be  present  on  either  side  of  a  plate.  In  the  example  in  Figure 
3-13,  a  total  of  17  nodes  is  assumed.  Nine  of  these  are  h^^  apart  in  the 
plate  1  region,  and  the  remaining  eight  are  h2  apart  in  the  plate  2  region. 
Nodes  1,  2,  8  and  9  are  'false'  nodes  corresponding  to  the  portion  of 
the  fastener  in  plate  1,  and.nodes  lO,  11,  16  and  17  are  'false'  nodes 
Corresponding  to  >  the  plate  2  region. 

In  a  symmetric  double  lap  joint . (see  Figure  3-12),  the  fastener  is 
discretized  only  up  to  the  plane  of  sy^etry  (midplane  of  plate  2).  In 
this  case,  m  in  Figure, 3-13  represents  half  of  the  number  of  plies  in 
plate  2. 

In  the  example  in  Figure  3-13,  m  +  n  +  8  nodes  are  present.  Of  these, 
m  +n  are  'physical'  nodes,  and  8  are  'faxse'  nodes.  The  enforcement  of 
the  governing  equation  at  the  ’physical'  nodes  yields  m  +  n  equations. 

The  imposition  of  the  boundary/continuity/symmetry  conditions  yields  8 
equations  (see  Section  3.6).  These  m  +  n  +  8  equations  provide  the 
solutions  for  the  m  +  n  +  8  nodal  displacements. 

Th."'  properties  of  a  ply  are  represented  by  an  equivalent  stiffness 
at  the  corresponding  node  in  this  approach.  The  computed  nodal  displace¬ 
ments  are  used  in  the  calculation  of  the  average  ply  loads.  When  a  ply 
failure  is  predieted,  the  ^.corresponding  nodal  stiffness  is  modified  to 
appropriately  represent  this  failure. 
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3.7.2  Central  Difference  Expressions 


0. 
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Assume  that  the  nodes  in  a  region  are  separated  by  a  distance  h.  In 
this  case  Az  ■  h  for  the  region,  and  the  derivatives  of  the  nodal  displace¬ 
ments  with  respect  to  z  may  be  approximated  by  the  following  central 
difference  expressions,  at  node  i: 


2hu^  ■ 

-u.  ,  +  u. 

i-1  i 

9  *» 

h  u .  - 

u.  -  -  2u 

i 

1-1 

Ilf 

2h  u^ 

“  ■  “l-2 

4  *  •  * » 

h 


i+1 


^“i-1  “  ^“i+1  “l+2 


“l-2  “i-1  ^“i  “l+l  “i+2 


(3-18) 


The  error  in  the  above  approximations  is  of  the  order  of  h  .  In  the  example 
in  Figure  3-13,  plate  1  is  divided  into  n  plies  that  are  h^^  in  thickness, 
and  plate  2  is  divided  into  m  plies  that  are  h2  in  thickness,  h^  and 

must  be  of  the  same  order  of  magnitude. 

It  is  noted  tliat  the  difference  approximations  of  derivatives  to  the 
fourth  order  (Equations  3-18)  involve  a  maximum  of  two  codes  on  either  side 
of  the  node  where  the  derivative  is  approximated.  The  enforcement  of  the 
governing  equation  at  every  physical  node,  therefore,  requires  two  'false* 
nodes  on  either  side  of  each  plate  (see  Figure  3-13).  The  boundary/continuity/ 
symmetry  conditions  are  enforced  at  the  first  or  the  last  nodal  location  in 
each  plate.  Referring  to  these  conditions  in  Equations  3-16  and  3-17,  and 
to  their  relationships  to  derivatives  to  the  third  order  of  the  nodal  displace¬ 
ment  (Equations  3-6  to  3-8),  it  is  seen  that  these  conditions  also  require  the 
two  'false'  nodes  on  either  side  of  each  plate.  A  general  node  arrangement 
for  an  n-ply  plate  bolted  to  an  m-ply  plate  is  shown  in  Figure  3-lA.  A  total 
of  8  'false*  nodes  are  required  for  the  discussed  finite  difference  formulation. 

The  'false*  nodes  are  assumed  *o  be  mirror  reflections  about  the  first 
physical  node  at  every  boundary.  Referring  to  Figure  3-14,  node  2  Is  assumed  to 
have  the  properties  (stiffness)  of  node  4,  and  node  1  Is  assumed  to  be  Iden- 
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tical  to  node  5.  Likewise,  nodes  nfS,  n+4,  n4-5,  n+6,  m+n+7  and  in+n+8  are 
mirror  reflections  of  nodes  n+1,  n,  n+9,  n+8,  m+n+5  and  m+n+4,  respec¬ 
tively. 


A  central  difference  approximation  of  the  governing  equation  is  obtained 
by  incorporating  Equation  3-18,  and  similar  difference  expressions  for 
Equation  3-11,  into  Equation  3-10.  The  foUowing  difference  expressions  are 
used  in  the  process: 


h^u^^  -  u^_2  -4  u^_j^  +  6u^  -4  u^^^  +  u^_^2 


qi  *  -k^u^  -  k^u^ 


(3-19) 
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‘*i-l  ‘2  qi  +  q 


i+1 


^■*^l-l“i-l“^l-l“l-l^"^^“Vi"Vi^  ^"\+l“i+l  "^1+1“ i+1^ 


“‘'l-l“i-l  ^*'l“l"'‘i+l"l+l  "  ^i-l“i-l  V  i~  ^i+l“l+l 


Substituting  the  above  expressions  into  Equation  3-10,  the  governing 


equation  at  the  i^^  node  becomes: 


“1-2 - * -!5p>  “i-i*<‘*T|r  *  — >  “i 


“1«  *  “1«  -  <4^ 


"i-1 


9  h„  2_  _  ^a^i+1  — 

-<  k  *  -ii’  "A"!  *  <  >"i« 


(3-20) 


a  -1  for  plate  1,  and  a  -  2  for  plate  2.  For  the  1*^^  ply,  the  values  of 
k,k  and  u  are  determined  by  its  state  of  damage  (see  Equation  3-11). 


Equations  3-6  to  3-8  express  the  fastener  bending  moment,  transverse 
shear  and  the  bending  slope,  respectively,  in  terms  of  u  and  q  derivatives. 
Central  difference  expressions  for  M,  y  and  yield  the  following  boundary/ 
continuity /symmetry  conditions : 


At  z»0  (fastener  head  location/  top  of  the  joint) 
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At  z*tj^  (joint  interface  between  plates  1  and  2) 
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At  z»t2^+t2  (fastener  nut  location/bottom  of  the  joint) 


(vii)  V2  ■  0: 
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(3-27) 


(vill)  M2  -  R2’P2 
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XGA  '  “ttfm+7  “  2"n+nri-8  2  n+mf6  n+m+6 


(3-28) 


3.8  Solution  Procedure 


The  goveri.xng  Equation  3-20  is  enforced  at  the  n+ci  'physical'  nodes. 
Equations  3-21  to  3-28  are  enforced  at  the  bounder,  nodes  (identified  in 
Figure  3-14  as  nodes  3  and  n+2  in  plate  1,  and  nodes  n+7  and  m+n+6  in 
plate  2).  This  provides  a  system  of  lu+r.+S  equations  that  determine  the 
through-the-thickness  fastener  displacement.  Figure  3-13  shows  a  matrix 
representation  of  the  problem  defined  in  Figure  3-13  (n=5,m=4).  The 
functions  f(u^)  are  set  equal  to  zero  until  the  ply  corresponding  to  mode 
i  is  damaged. 

The  system  of  in+n+8  equations,  in  general,  is  solved  for  the  nodal 
displacements  by  using  a  standard  matrix  decomposition  or  matrix  inversion 
routine.  It  is  seen,  from  Figure  3-15,  that  the  coefficient  matrix  to  be 
Inverted  is  banded  in  nature.  A  special  purpose  Gaussian  triangularization 
computer  code  was  developed  to  solve  for  the  nodal  displacements.  When 
the  coefficient  matrix  size  is  large,  this  provides  an  economical  means  of 
obtaining  solutions,  especially  when  the  solution  procedure  is  repeated 
many  times  to  predict  progressive  failures  in  the  bolted  plates. 

The  finite  difference  formulation  of  the  fastener  analysis  has  been 
programmed  to  be  the  FDFA  co.jputer  code.  Convergence  studies  on  the  analysis, 
correlation  of  FDFA  predictions  with  available  analytical  solutions,  and 
sample  predictions  using  FDFA  are  presented  below. 

3.9  Convergence  Study  on  the  Fas : ener  Analysis  (FDFA) 

The  number  of  actual  plies  in  a  bolted  laminate  determines  the  number  of 
nodes  in  the  portion  of  the  fastener  within  that  laminate.  A  physical  ply  can 
also  be  divided  equally  into  two  or  more  plies  of  smaller  thicknesses,  to 
improve  the  accuracy  of  the  solution.  Referring  to  Figure  3-13,  If  plate  1  is 
a  laminate  with  n  plies,  the  portion  of  the  fastener  within  plate  1  can  be 
divided  into  £sn  segments  that  are  t^/(£zn)  in  thickness,  where  I  is  any 
integer  ^1.  In  Figure  3-13,  n  ■  5  and  £  ■  1. 

The  effect  of  the  number  of  nodes  per  bolted  plate  on  the  displacement 
solution  was  studied  by  considering  a  ateel-to-steel  bolted  Joint  example. 

The  steel  plates  were  0.125  inch  in  thickness,  and  transferred  a  tensile  load 


Figure  3-15.  The  Matrix  Equations  for  n=5  and  m=4  (see  Figure  3-13). 


In  a  single  lap  configuration  via  a  0.25  inch  diameter  steel  bolt.  The  steel 
bolt  was  assumed  to  have  a  modulus  (E)  of  30  x  10^  and  a  Poisson's  ratio  (2) 
of  0.30.  The  steel  plates  were  assumed  to  have  an  effective  modulus  of 
7.6S  X  10^  psi,  and  the  bolt  was  assumed  to  be  free  to  rotate  at  the  head  and 
nut  locations  (zero  moment  or  pin  conditions).  Figure  3-16  shows  the 
predicted  fastener/plate  displacement  variation  through  the  thickness  of 
plate  1,  for  an  applied  load  of  1000  lbs,  when  the  number  of  nodes  in  plate  1 
Is  increased  from  5  to  40.  It  Is  seen  that  beyond  20  nodes,  the  displacement 
distribution  is  relatively  unaffected. 

Figure  3-17  presents  the  fastener /plate  displacement  distribution  through 
the  thickness  of  plate  1,  when  the  1  kip  load  transfer  is  affected  in  a 
double  lap  configuration.  In  this  case,  the  properties  were  assumed  to  be 
the  same  as  those  in  Figure  3-16,  and  plata  2  was  0.25  inch  In  thickness.  The 
results  in  Figure  3-17  also  indicate  that  beyond  20  nodes,  the  displacement 
solution  is  relatively  unaffected. 

It  is  therefore  necessary  to  divide  a  metallic  plate  Into  at  least  20 
identical  "plies"  to  obtain  a  converged  displacement  solution,  which  ensures 
an  accurate  prediction  of  the  load  distribution  in  the  thickness  direction. 

In  laminated  (non-metallic)  plates,  a  minimum  of  30  nodes  per  plate  is 
required  to  obtain  accurate  solutions. 

3. 10  Comparison  of  FDFA  Predictions  with  Available  Analytical  Solutions 

Harris  and  OJalvo  (Reference  3-1)  have  obtained  analytical  solutions  for 
metal-to-metal  load  transfer  in  single  and  double  shear  configurations,  Khen 
the  bo\indary  conditions  at  the  fastener  head  and  nut  Ijcations  are  restricted 
to  zero  moment  (free  pin)  or  zero  slope  (rigid  head  or  lut).  Predictions 
made  by  the  finite  difference  fastener  analysis  (FDFA)  ode  in  SASCJ  were 
compared  with  the  results  in  Reference  3-1  to  establish  the  validity  of  FDFA. 

The  first  example  selected  for  the  correlation  study  is  a  single  shear 
situation  with  no  rotational  constraint  (zero  moment  or  free  pin  condition) 
at  the  fastener  head  and  nut  locations.  The  bolted  plates  are  0.125  inch 
thick  steel  plates,  and  are  assumed  to  have  an  effective  modulus  of  7.68  x  10^ 
psi.  The  steel  fastener  is  0.25  inch  in  diameter,  with  E  ■  30  x  10®  psi 
and  V  ■  0.3.  The  fastener  is  assxnned  to  be  a  flexural  beam  in  one  case  (rigid 
in  shear),  a  shear  beam  in  the  second  case  (rigid  in  bending),  and  a  completely 
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Displacement  Solution  for  a  Single  Shear  Situation. 


Double  Shear  Configuration 


rigid  Dean  in  the  third  case.  For  the  flexural  beam,  the  analysis  in  Reference 

3-1  set  GA  to  be  infinite  in  magnitude,  and  FDFA  assumed  XgA  to  be  10^^  lbs. 

For  the  shear  beam,  FDFA  assumed  El  (infinity  in  Reference  3-1)  to  be 
10  2 

10  Ib-ln  .  For  the  completely  rigid  beam,  the  analysis  in  Reference  3-1 

12 

set’El  and  GA  to  be  infinity,:  and  FDFA  set  E  to  be  3  x  10  psi.  A  comparison 
between  the  predictions  made  by  FDFA  and  the  analysis  in  Reference  3-1  is 
presented  in  Table  3-1.  Refer  to  Figure  3-11  for  the  definition  of  z/t^,  where 
t^  is  the  thickness  of  the  top  plate  (plate  1),  and  z  is  measured  from  the 
fastener  head  location  (top  of  plate  1)  toward  the  nut  location.  For  the 
steel  plates,  the  loading  rate  (q  in  lbs /inch)  is  expressed  by  q  >-ku 
(see  Equation  3-11),  where  k  is  the  effective  modulus  of  the  plate  material 
and  u  is  the  fastener/plate  displacement.  The  quantity  qt^/P  represents  a 
load  interaction  parameter  which  piroyides  a  measure  of  the  deviation  in  the 
through-the-thickness  load  distribution  from  the  uniform  value  of  unity 
corresponding  to  a  two-dimensional  analysis  of  the  plate.  P  is  the  applied 
load  on  the  steel  plates  (see  Figure  3-11). 

Table  3-1  indicates  that  the  FDFA  predictions  are  within  5  percent  of  the 
analytical  predictions  made  In, -.Reference  3-1.  It  is  noted  that  rigidity  in 
shear  or  bending,or  both, is  only  approxinated  by  the  FDFA  code.  The 
represents ''ion  of  the  rigidities  by  numbers  larger  than  those  shown  in 
Table  3-1  result  in  numerical  problems  that  result  from  an  ill-conditioned 
coefficient  u> ''>r  the  unknown  displacements.  Nevertheless,  FDFA  predic¬ 
tions  and  those  bas^  '  the  analysis  in  Reference  3-1  demonstrate  an 
excellent  agreement  (erroi  5Z),  validating  the  FDFA  code.  It  is  also  noted 
that  the  three  types  of  rigidity  yield  approximately  the  same  results  using 
either  analysis. 

sucuuu  example  selected  for  the  correlation  study  is  the  same  as  the 
first  example,  but  the  fastener  head  and  nut  locations  are  assumed  to  be 
rigidly  constrained  (zero  bending  slope).  Table  3-2  presents  the  comparison 
between  predictions  made  by  FDFA  and  the  analysis  in  Reference  3-1.  Again, 
the  difference  between  the  two  sets  of  results  is  less  than  5  percent  and  the 
two  types  of  rigidities  yield  approximately  the  same  results. 

The  last  example  considered  for  the  validation  of  the  FDFA  code  is  a 
double  shear  situation  where  two  0.0625  inch  thick  steel  plates  (labeled  1 
in  Figure  3-12)  transfer  the  applied  load  to  a  0.125  inch  thick  steel  plate 
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Steel  plates  1  and  2  are  0.125  Inch  in  thickness,  with  an  effective  modulus  of  7  '8  x  1*'^  psi. 
The  steel  fastener  is  0.25  inch  in  diameter,  with  E  =  30  x  10^  psi  and  v  =  0.3. 


(labeled  2  in  Figure  3-12)  via  a  0.25  inch  diameter  steel  fastener.  The 

fastener  is  assumed  to  be  completely  rigid  (El  and  GA  •*•“),  and  is  assumed  to  be 

free  to  rotate  at  the  head  and  nut  locations  (zero  moment  conditions).  The 

FDFA  code  represented  a  rigid  fastener  by  increasing  its  modulus  from 
6  10 

30  X  10  to  30  X  10  psi.  The  analysis  in  Reference  3-1  predicted  a  uniform 
load  distribution  (qt^/P  =  1  for  0  <  zlt^  <  1).  The  FDFA  code  also  predicted 
a  uniform  load  distribution  (qt^/P  =  1.034  for  0  <  z/t^  <  1),  and  indicated 
a  3%  error  in  the  recovery  of  the  applied  load. 

The  FDFA  code,  therefore,  predicts  accurately  the  closed  form  solutions 
for  the  discussed  joint  situations.  The  error  in  the  FDFA  predlctlon.s  is 
restricted  to  within  5  percent,  and  the  FDFA  code  is  also  applicable  in 
situations  where  closed  form  solutions  are  difficult  to  obtain.  These 
Include  non-metallic  bolted  plates,  and  general  elastic  constraints  at  the 
fastener  head  and  nut  locations. 

3. 11  Effect  of  Single  and  Double  Shear  Loan  Transfer  on  the  Through-the-Thlckness 
Load  Distribution 

The  difference  in  the  through-the-thlckness  load  distributions  between 

single  and  double  shear  load  transfer  situations  is  demonstrated  in  Figure  3-18. 

In  this  case,  the  bolted  plates  are  made  of  steel,  and  the  1/4  inch  diameter 

steel  fastener  is  assumed  to  be  free  to  rotate  at  the  head  and  nut  locations. 

* 

The  load  Intensity  in  a  ply  Is  quantified  by  -ku  (Ibs/in),  and  the  load 
intensity  for  a  uniform  distribution  is  expressed  as  (Ibs/ln)  for  the  top 

plate  (plate  1).  The  load  Intensity  ratio  -  kut^^/P,  therefore,  is  a  measure 
of  the  variation  from  the  uniform  distribution.  Figure  3-18  indicates  that  the 
load  intensity  ratio  varies  between  -2  and  4  for  the  single  shear  situation, 
and  Is  approximately  1  (near  uniform  distribution)  for  the  double  shear 
situation. 

If  the  fastener  head  and  nut  locationr  are  constrained  against  rotation 
(zero  bending  slope),  the  single  shear  situation  also  results  in  a  near  uniform 
load  distribution  (see  Table  3-2).  In  practical  joints,  rotational  constraints 
at  the  fastener  head  and  nut  locations  are  introduced  by  the  washers,  and  are 
influenced  by  the  fastener  torque  (aee  Section  3.6).  Regardless  of  what  the 
empirical  elastic  constraint  equations  are  at  these  locations  (M  ■  Ri^),  the 
load  intensity  ratio  will  be  approximately  unity  in  practical  situations,  and 


63 


will  not  vary  as  much  as  in  the  situation  considered  in  Figure  3-18.  An 
exception  is  a  flush-head  or  countersunk  fastener  in  which  the  head  is 
practically  unconstrained  against  rotation.  But,  the  FDFA  code  is  currently 
restricted  to  protruding  head  fasteners. 

3.12  Effect  of  Fastener  Size 

The  effect  of  fastener  size  on  the  through-the-thlckness  load  distribution 
is  shown  in  Figure  3-19.  The  example  assumes  two  0.125  inch  thick  steel  plates 
to  be  joined  by  a  steel  fastener  in  a  single  shear  configuration.  The  head  and 
nut  locations  of  the  fastener  are  assumed  to  be  unrestrained  (zero  momenr 
conditions).  Figure  3-19  indicates  that  for  D/tj^  values  that  are  greater  than 
or  equ£d  to  unity,  the  fastener  size  has  a  negligible  effect  on  the  through- 
the-thlckness  load  distribution.  For  D/tj^  _>  1,  the  load  intensity  ratio  varies 
from  -2  to  4  as  z/t^  varies  from  0  to  1.  But,  for  D/t^  values  that  are  less 
than  unity,  the  fastener  bending  and  shear  effects  Increase  with  a  decrease  in 
D/tj^*  For  example,  the  load  intensity  ratio  varies  from  near 

0  to  10  as  z/t^  varies  from  0  to  1.  In  these  situations,  the  large  fastener 
displacements  generally  result  in  fastener  failures. 

3.13  Effect  of  Fastener  Material 

Three  fastener  materials  were  considered  for  this  study  —  steel, 
titanltm  and  aluminum.  The  effect  of  the  material  on  the  through-the-thickness 
load  distribution  (load  intensity  ratio)  is  presented  in  Figures  3-20  to  3-22. 
When  B/t^  is  greater  than  or  equal  to  unity,  the  material  has  a  negligible 
effect  on  the  load  distribution.  In  this  case,  the  fastener  behaves  like  a 
rigid  fastener,  and  the  differences  among  the  moduli  of  the  three  uatcrials 
have  no  effect  on  the  load  distribution.  But,  when  D/t^  is  less  than  unity 
(1/2  in  Figure  3-21,  and  1/4  in  Figure  3-22),  the  fastener  flexlbllitv  is 
Influenced  significantly  by  the  modulus  of  the  material.  In  this  case,  the 
material  with  the  lowest  modulus  (alimlnum  in  the  chosen  example)  yields  the 
largest  load  intensity  ratio. 

3*1^  Effect  of  Stacking  Sequence  on  Throuch-the-Thickness  Load  Distribution 

When  one  of  the  bolted  plates  is  a  laminate,  the  stacking  sequence  in¬ 
fluences  the  ply  load  distribution  and  consequently  the  joint  strength.  To 
Illustrate  this,  selected  compos it e-to-alumlnum  tests  from  Reference  3-6  arc 
considered.  The  selected  test  cases  Involve  50/40/10,  ASl/3501-6  graphite/ 
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epoxy  laminates  of  three  different  stacking  sequences.  The  laminates  are 
0.119  inch  in  thickness,  contain  20  plies,  and  are  bolted  to  0.31  inch  thick 
aluminum  plates  with  0.25  inch  diameter,  protruding  head,  steel  fasteners. 

The  load  transfer  occurs  in  a  single  shear  configuration  (Figure  3-11). 

The  variation  in  the  ply  loads  with  a  change  in  the  laminate  stacking 
sequence  was  predicted  using  the  FDFA  code.  The  aluminum  plate  was  divided 
into  30  layers  of  identical  pruperties  (E  ■  10  Msi,  v  ■  0.3) .  Each  ply  in 
the  laminate  was  divided  equally  into  two  layers  of  identical  properties.  The 
20T>ply  laminate  was,  i  nerefore,  modeled  as  a  40-ply  plate  with  the  following 
properties:  E^^j^  ■  18.3  Me i,  Msi,  Vj^2  “  Msi.  The 

aliminum  and  composite  plates  were  assumed  to  be  2  inches  long  and  1  inch  wide, 
with  the  0^25  inch  diameter  hole  at  the  center  (E/D  ■  W/D  ■  4).  The  effective 
foundation  modulus  for  each  ply  type  (fiber  orientation)  was  computed  as 
■described  in  Section  3.5.  This  yielded  an  effective  foundation  modulus  (k)  of 
28.8  Msi  for  alum  intun,  30.5  Msi  for  the  0**  ASl/3501-6  ply,  19.2  Msi  for  the 
+  or  -45°  ASl/3501-5  ply,  and  ^.76  Msi  for  the  90°  ASl/3501-6  ply.  The 
elastic  constraints  at  the  fastener  head  and  nut  locations  were  assumed  to  be 
of  the  form  Shown  in  Equation  3-16,  and  the  disk  radius  in  Equation  3-15  was 
assumed  to  be  1/8  inch.  ^^^  ^  ^  ^  ^ 

Using  the  above  input  data,  FDFA  predicted  the  ply  load  distributions  in 
50/40/10  layups  of  three  different  stacking  sequences  —  [(4 5/0/ -4 5/ 0)2/0/90]^, 
[02/+45/02/445/90]^  and  [90/+452/Oj3|j*  Figure 3-23  shows  how  different  the 
ply  loads  are  when  the  stacking  sequence  is  varied.  This  also  has  a  direct 
influence  on  the  Joint  strength. 
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Figure  ■’-23.  Ply  Loads  In  ASl/3501-6  Graphite/Epoxy  Laminates  with  a  50/40/10 
Layup  and  Different  Stacking  Sequences,  When  Bolted  to  Aluminum 
Plates  in  a  Single  Shear  Configuration. 
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Figure  3-23.  Ply  Loads  in  AS1/3S01-6  Craphlte/Epoxy  Laminates  with  a  50/40/10 
Layup  and  Different  Stacking  Sequences,  When  Bolted  to  Aluminum 
Plates  In  a  Single  Shear  Configuration  (Continued). 


Figure  3-23,  Ply  Loads  in  ASl/3501-6  Graphite/ Epoxy  Laminates  with  a  50/40/10 
T>ayup  and  Different  Stacking  Sequences,  When  Bolted  to  Aluminum 

Plates  In  a  Single  Shear  Configuration  (Concluded). 


SECTION  4 

STRENGTH  ANALYSIS  OF  SINGLE  FASTENER  JOINTS  IN  COMPOSITE  STRUCTURES 
4.1  Introduction 

The  two-dimensional  anisotropic  plate  analysis  described  in  Section  2 
(FIGEOM)  and  the  finite  difference  fastener  analysis  described  in  Section  3 
(FDFA)  are  incorporated  into  a  failure  procedure  to  develop  a  strength 
analysis  for  single  fastener  joints  in  composite  structures.  A  general 
fastener  location  in  a  bolted  structure  (see  Figure  1-1)  Involves  the  type 
of  loading  shown  in  Figure  4-1.  The  general  situation  can  be  analyzed  as  a 
superposition  of  an  unloaded  hole  situation  and  a  fully-loaded  hole  situation, 
as  shown  in  Figure  4-1.  The  unloaded  hole  case  is  analyzed  using  the  two- 
dimensional  plate  analysis  (FIGEOM),  and  does  not  involve  the  fastener  analysis 
(FDFA).  The  fully-loaded  hole  situation  Involves  both  the  analyses,  and 
Section  4-2  describes  the  failure  procedure  for  this  case. 

The  loaded rhole  si cuatlon  is  analyzed  using  a  progressive  failure 
procedure  that  predicts  local  ply  failures  and  delamlnations  until  the  bolted 
plate  cannot  carry  any  additional  applied  load.  The  employed  ply  failure 
criteria  and  the  delamlnatlon  criterion  are  discussed  in  Sections  4.4  and 
4.5,  respectively. '  The  strength  analysis  has  been  programmed  to  be  the 
SASCJ  (Strength  Analysis  of  Single  Fastener  Composite  Joints)  computer  code. 
Sample  strength  predictions  using  SASCJ  are  presented  in  Section  5. 

4*  2  Strength  Analysis  Procedure'  for  Fully-Loaded  Holes 

The  strength  of  laminates  with  fully-loaded  holes  is  predicted  using 
the  procedure  outlined  in  Figure  4-2.  A  two-dimensional  stress  analysis 
(FIGEOM),  accounting  for  finite  dimensions  of  the  bolted  plates,  is  initially 
performed  on  each  bolted  plate.  Computed  plate  stresses  are  used  to  calculate 
the  effective  moduli  pf  the  various  ply  types  in  each  bolted  plate  (see 
Section  3.5).  The  inplane  strains  computed  by  the  FIGEOM  code  are  used  to 
obtain  the  stress  state  in  each  ply.  The  ply  stresses  around  the  hole  boundary 
are  integrated  to  yield  the  bearing  load  in  each  ply  (see  Section  3.5).  The 
Inplane  stresses  in  each  ply,  per  unit  bearing  load,  are  incorporated  into 
selected  failure  criteria  to  compute  the  ply  (bearii^)  loads  corresponding 
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Figure  4-1.  Schematic  Representation  of  a  General  Single  Fastener 
Situation  as  a  Superposition  of  Unloaded  and  Fully- 
Loaded  Hole  Situations. 
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Flowchart  for  the  Strength  Analysis  of 
Laminates  with  Fully-Loaded  Holes. 


to  the  various  inplane  failure  modes.  The  selected  failure  criteria  are  pre¬ 
sented  in  Section  4.4. 

The  effective  moduli  and  the  ply  bearing  loads  corresponding  to  the 
various  failure  modes,  for  all  the  plies  in  each  bolted  plate,  are  incor¬ 
porated  into  the  fastener  analysis  described  in  Section  3,  The  Initial 
fastener  analysis  on  the  undamaged  plates  assumes  the  u  terms  for  all  the 
plies  to  be  zero  (see  Equations  3-"'“  to  3-28),  and  computes  the  distribu¬ 
tion  of  the  applied  bearing  Ir  •.  among  the  various  plies.  Comparing  these 
'■!”  loads  with  the  stored  .illure  values  for  inplane  ply  failures,  the 
joint  loau  _  -■'•--n  nuing  to  the  earliest  ply  failure  is  obtained.  The 
fastener  analysis  also  computes  approximate  shear  strain  values  at  the  inter- 
facial  locations  between  adjacent  plies.  Incorporating  these  into  an  inter¬ 
laminar  failure  criterion,  the  joint  load  corresponding  to  the  earliest  inter¬ 
laminar  failure  (delamination)  is  obtained.  The  smaller  of  the  two  joint  loads, 
corresponding  to  the  earliest  Inplane  and  Interlaminar  failures,  determines 
the  first  fall.ire  in  a  bolted  plate  and  the  corresponding  joint  load  value. 

The  effective  moduli  of  the  damaged  plies  are  reset  to  appropriately 
represent  the  predicted  failure  modes.  The  revised  moduli  are  Incorporated 
into  the  fastener  analysis,  and  the  procedure  is  repeated  to  predict  the 
next  failure  mode  and  the  corresponding  Joint  load.  When  any  ply  is  pre¬ 
dicted  to  fall  totally,  the  analysis  computes  the  redistribution  of  the 
corresponding  joint  load  among  the  remaining  effective  plies,  and  determines 
if  any  other  concomitant  ply  failure  is  precipitated,  This  process  is 
repeated  until  one  of  the  bolted  plates  becomes  ineffective  in  transferring 
the  applied  load  (joint  failure). 

The  SASCJ  computer  code  is  restricted  to  protruding  head  type  fasteners, 
and  assumes  that  fastener  failure  is  precluded.  It  can  analyze  any  combination 
of  laminated  and  metallic  plates,  bolted  together  in  a  single  lap  or  double 
lap  configuration. 

4 . 3  Strength  Analysis  Procedure  for  Partially-Loaded  Holes 

A  general  fastener  location  in  a  bolted  plate  transfers  a  fraction 
(a)  of  the  total  applied  load  via  the  fastener,  the  remainder  (1-a)  being 
by-passed  to  the  next  fastener  location  (see  Figures  1-1  and  4-1).  In  this 
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casa,  the  stress  state  at  the  fastener  location  is  computed  as  a  superposition 
of  the  stress  states  corresponding  to  the  unloaded  and  fully-loaded  hole 
situations.  Figure  4-3  presents  a  schematic  representation  of  how  the  averaged 
stresses  are  obtained  to  predict  net  section,  shear-cut  and  bearing  failures 
in  the  plies  using  average  stress  failure  criteria  (see  Section  4.4).  For  a 
unit  applied  load,  the  averaged  stresses  in  the  laminate  with  an  unloaded 
hole,  when  subjected  to  a  load  of  (1-a),  and  the  averaged  stresses  in  the 
laminate  with  a  fully  loaded  hole,  when  subjected  to  a  load  of  a,  are  com¬ 
puted  separately  and  added.  Incorporating  the  combined  averaged  stresses 
into  the  appropriate  failure  criteria,  the  applied  load  corresponding  to  a 
ply  failure  is  computed. 

In  the  case  of  fully-loaded  holes,  progressive  failure  prediction 
Involves  the  repetition  of  the  fastener  analysis  with  revised  ply  properties 
after  every  ply  failure.  The  two-dimensional  analysis  (FIGECM)  is  only 
carried  out  once  (see  Section  4.2).  But,  in  the  case  of  partially-loaded 
holes,  a  ply  failure  will  affect  the  unloaded  and  the  fully-loaded  hole  con¬ 
tributions  to  the  local  stresses.  Hence,  progressive  failure  prediction  in 
the  partially-loaded  case  involves  repeating  FIGEOM  and  FDFA  analyses  after 
total  ply  failures. 

4,4  Inplane  Failure  Criteria  ■ 

The  SASCJ  code  permits  the  user  to  select  any  of  the  following  five 
failure  criteria  for  the  prediction  of  ply  failures  based  on  Inplane  stresses 
and  strains:  (1)  point  stress  failure  criterion,  (2)  average  stress  failure 
criterion,  (3).  maximum  (fiber  directional)  strain  criterion,  (4)  Hoffman 
criterion,  and  (5)  Tsal-Hillcriterlbn,  The  first  two  criteria  predict 
three  modes  of  failure  in  each  ply— net  section,  shear-out  and  bearing.  The 
maximum  strain  criterion  predicts  ply  failure  based  on  fiber  directional 
strain.  The  Hoffman  and  Tsai-Hlll  criteria  predict  ply  failure  accounting 
for  biaxial  stress  interaction  that  is  ignored  by  the  first  three  criteria. 

The  point  stress  failure  criterion  predicts  net  section,  shear-out  and 
bearing  failures  when  the  appropriate  stress  components  at  selected  locations 
attain  unnotched  specimen  failure  values  (see  Figure  4-4;.  a°^,  a^°  and  a^^^ 

■Xv...  .  O  .  O'  O 

are  called  characteristic  distances  (see  Reference  4-1).  When  0^(0,  D  +a°®) 
exceeds  the  unnotched  tensile  dr  campresslve  strength  of  the  ply,  as  appropriate. 


The  value  at  this  location  determines 
net  section  failure 


■The  T  value  at  this 
xy 

location  determines  shear- 
out  failure 


The  a  value  at  this 


location  determines 
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Figure  4-4.  The  Characteristic  Distances  used  in  the  Point 
Stress  Failure  Criteria 


a  net  section  ply  failure  is  predicted.  When  a  (D  +  a^^®,  0)  exceeds  the  un- 

3C  O 

notched  compressive  strength  of  the  ply,  a  bearing  mode  of  ply  failure  is 
predicted.  When  (a^  ,  D/2)  exceeds  the  unnotched  ply  shear  strength,  a 
shear-out  mode  of  ply  failure  is  predicted.  The  average  stress  failure 
criterion  predicts  these  failures  based  on  averaged  values  of  the  mentioned 
stress  components  over  selected  characteristic  distances  (d'^^,  d^°,  and  d^^®) 
that  are  larger  in  magnitude  compared  to  those  used  in  conjunction  with  the 
point  stress  criterion  (see  References  4-2  and  4-3,  and  Figure  4-5) 

Of  the  three  ply  failure  modes,  only  the  net  section  mode  causes  the 
ply  to  become  almost  Ineffective  (total  failure).  The  bearing  mode  of 
failure  causes  the  ply  to  suffer  a  reduction  in  its  effective  modulus 
without  losing  its  load-carrying  capacity.  The  shear-out  mode  of  failure 
causes  a  ply  to  become  Ineffective  only  when  it  Is  delaminated  from  the 
adjacent  plies  (see  Section  4.5).  When  a  ply  suffers  any  of  the  above 
failures,  its  load  versus  deflection  response  Is  at  the  knee  of  the  bilinear 
representation  In  Figure  3-8.  The  damaged  ply  can  carry  additional  load 
until  total  ply  failure  is  precipitated.  The  SASCJ  computer  code  auto¬ 
matically  stores  the  damage  state  in  every  ply  in  the  bolted  plates,  and  re¬ 
assigns  values  for  ply  moduli  to  appropriately  represent  predicted  ply 
failures.  When  a  ply  suffers  total  failure,  its  modulus  is  set  equal  to 
zero  (see  Equation  3-11),  and  the  redistribution  of  the  joint  load  among 
the  remaining  plies  is  computed.  A  typical  overall  load  versus  deflection 
behavior  of  the  joint  is  sho\.m  in  Figure  4-6,  indicating  the  effects  of  local 
and  total  ply  failures. 

The  maximum  strain  (fiber  directional),  Hoffman  and  Tsal-Hlll  criteria 
are  applied  along  a  path  that  is  concentric  to  the  fastener  hole,  at  a 
characteristic  di.«»tance  (a^)  from  the  hole  boundary  (see  Figure  4-7  and 
Reference  4-4),  The  location  along  this  path  where  the  selected  criterion 
is  satisfied  determines  the  failure  location..  The  maximum  strain  criterion 
predicts  fiber  failure  in  a  ply  (total  ply  failure)  when  its  fiber  directional 
strain  exceeds  the  failure  vaiue  (e^^  or 
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The  Hoffman  failure  criterion,  based  on  inplane  ply  stresses,  states 
that  total  ply  failure  will  occur  when  the  failure  index  (H)  In  the 
following  equation  reaches  a  value  of  unity: 

. . .  ..  2 .  .  .  . 


2  2 

O^/S^  -  H 


(i-1) 


In  the  above  equation,  0^,  and  are  the  ply  stresses  in  the  fiber 
coordinate  system,  and  are  the  uniaxial  tensile  and  compressive 
material  strengths  along  the  fiber  direction  (1),  and  Y  are  the  uniaxial 
tensile  and  compressive  material  strengt’  ,<  perpendicular  to  the  fiber 
direction  (2),  and  S  is  the  material  shear  strength  in  the  1-2  plane. 


In  the  SASCJ  code,  the  Hoffman  criterion  is  applied  along  a  path  that 
is  concentric  to  the  fastener  hole,  defined  by  the  characteristic  distance 
a^  (see  Figure  4-7).  At  selected  points  along  this  path,  the  following 
expressions  for  the  failure  value  of  the  ply  load  (?£)  are  computed: 


The  location  where  the  smallest  non-negative  value  fox  is  computed 
identifies  the  failure  initiation  point. 

The  Hoffman  criterion  predicts  total  ply  failure  and  the  failure 
location,  but  does  not  Identify  the  mode  of  failure.  The  failure  location, 
though,  generally  Indicates  the  possible  failure  mode.  Referring  to 
Figure  4-7,  if  failure  is  predicted  near  6“0°,  a  bearing  mode  of  failure 
is  suspected.  If  the  failure  location  is  near  0*90°,  a  net  section  mode 
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of  failure  is  suspected.  And,  Intermediate  values  of  6  indicate  a 
shear-out  mode  of  failure.  The  Tsai-Hill  criterion  can  be  obtained  from 
Equation  4-1  by  setting  and  “  Y^.  This  criterion,  therefore, 

does  not  account  for  different  strengths  under  tension  and  compression.  The 
ply  failure  load  (P^)  in  this  case  is  computed  to  be  1/^3”  (see  Equation  4-2 

A.''  Interlaminar  Failure  Criterion 

Delamination  between  plies  is  predicted  by  incorporating  computed 
shear  strains  at  the  interfaCial  locations  into  a  maximum  shear  strain 
criterion.  At  the  interface  betwen  plies  i  and  j,  for  example,  the  shear 
strain  is  computed  to  be: 

«  (u.  -  u.)  /h  (4-3) 

x2  I  j  a 

where  h^  is  the  ply  thickness  in  the  plate  containing  plies  i  and  j.  Thia 
expression  for  the  shear  strain  is  approximate.  Plies  1  and  j  are  assumed 
to  delaminate  when  exceeds  a  failure  value.  The  failure  value  for 

Yjjj,  will  be  determined  by  correlating  predictions  with  observations  for  a 
sample  test  case. 


SECTION  5 

STRENGTH  PREDICTIONS  USING  THE  SASCJ  CODE 
5.  1  Introduc  t Ion 

Tlie  developed  strength  analysis  for  laminated  and  metallic  plates 
(SASCJ  code)  requires  as  input  a  few  failure  parameters  and  post-failure  ply 
property  degradations  (see  Section  4).  Section  5.3  describes  how  these 
quantities  are  determined.  Subsequently,  the  validity  of  the  SASCJ  code  is 
established  by  comparing  its  strength  predictions  with  experimental  measure¬ 
ments  (see  Sections  5.4  to  5.6).  Test  results  from  Reference  5-1  are  used 
in  the  determination  of  the  mentioned  input  properties  and  for  establishing 
the  validity  of  the  SASCJ  code.  Figures  5-1  and  5-2  present  the  geometries 
of  the  bolted  laminated  and  t<etalllc  plates  used  in  Reference  5-1.  Table  5-1 
lists  the  various  static  tests  conducted  in  Reference  5-1. 

5 . 2  Initial  Studies  On  Failure  Parameters  For  Laminates  With 

Unloaded  (Open)  Holes 

The  strength  of  laminates  with  unloaded  or  loaded  holes  has  hitherto 
been  predicted  based  on  a  one-step  procedure.  The  laminate  or  the  individual 
lamina  stresses  or  strains  are  Incorporated  into  a  selected  failure  criterion 
to  predict  the  laminate  strength.  The  procedure  does  not  go  beyond  the  first 
ply  failure  when  lamina  stresses  or  strains  are  used  to  predict  laminate 
strength.  This  strength  prediction  procedure  results  in  failure  parameters  that 
are  dependent  on  the  laminate  layup,  and  requires  modlflcaiton  of  lamina 
properties  to  preclude  "non-critlcal"  failures.  Examples  are  presented  below. 

Consider  a  laminate  with  an  unloaded  (open)  5/16  in.  diameter  hole 
subjected  to  static  tensile  or  compressive  loading.  At  the  laminate  level, 
the  two-dimensional  laminate  stress  state  around  the  fastener  hole  is  in¬ 
corporated  into  a  selected  failure  criterion  to  predict  the  laminate  strength. 
The  point  stress  and  average  stress  failure  criteria  have  been  used  before  for 
this  problem  (see  References  5-2  to  5-4).  The  characteristic  distances  for 
net  section,  bearing  and  shear-out  failures,  respectively,  that  make  SASCJ 
predictions  agree  with  the  measured  tensile  strengths,  for  three  laminate 
layups,  are  presented  in  Figure  5-3.  'the  axial  tensile  and  compressive 
strengths,  and  the  shear  strengths,  of  the  laminate  layups  were  computed 
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Figure  5-2.  Metal  Plate  Dimensions  Used  in  Referenc 


TAPLE  5-1 .  STATIG  TESTS  CONDUCTED  IN  REFERENCE  5- 1 . 
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lARI.r.  5-1.  STATIC  VKSIS  CON'MTTED  IN  REFERENCE  5-1  (Continued) 
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Figure  5-3.  Failure  Parameters  for  the  Laminate  Level  Prediction  of 

Static  Tensile  Strengths  of  Laminates  with  Unloaded  (Open) 
5/16  Inch  Diameter  Holes. 
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usin,’,  lajnlnate  thcorV  :a  ASl/3501-f)  proporties  from  Keforence  5-5.  The 
cairactcristic  distances  that  make  SASCT  predictions  agree  with  the  measured 
C'-pressive  strengths  are  presented  in  Fieurc  5-4,  Predictions  in  both  cases 

,jre  sensitive  onlv  to  the  a‘  value,  since  the  laminates  fail  in  a  net 

■’■■pns  ■ 

section  mode.  The  a  value  for  tlie  average  stress  criterion  is  approximately 
2.5  times  the  value  for  the  point  stress  criterion,  as  observed  in  Reference  5-2. 
Note  tliat  a^^^^  for  the  highly  fiber-dominated  70/20/10  layup  is  much  larger  than 
that  Icr  the  50/40/10  and.  30/60/ 10  layups  under  tension. 

The  tensile  and  compressive  strengths  of  the  laminates  with  unloaded 
(open),  5/16  inch  di;imeter  holes  can  a:lso  be  predicted  at  the  lamina  level, 
based  on  first  ply  failurt%  Tlie  point  or  the  average  stress  criterion  (see 
Figures  t-H  and  -t-S),  the  Hoffman,  Tsai-Hill  or  the  maximum  strain  criterion 
(see  Figure  4-7)  may  be  used  at  the  ply  level.  Lami.nate  failure  is  assumed 
to  bo  precipitated  by  the  first  major  ply  failure.  When  the  Tsai-Hill  cri¬ 
terion  is  used,  matrix  failures  between  fibers  are  precluded  by  setting  the 

■  ■  4 

transverse  tensile  and  compre.ssive  strength  for  the  ply  to  be  large  (10  psi). 
When  the  maximum  strain  criterion  is  used,  only  the  fiber  directional  strain  is 
interrogated.  When  the  point  and  average  stress  criteria  are  used,  the 
strengths  of  the  +45°  plies  are  incrfeasCd  to  twice  the  actual  values,  and  the 
transverse  tensile  strength  is  increased  to' preclude  matrix  failures  between 
fibers  in  the  90°  plies;  Figures  5-5  and  5-6  present  the  characteristic  dis¬ 
tances  that  make  SASCJ  predictions  agree  accurately  with  the  measured  tensile 
and  compressive  strengths  of  three  laminates  with  unloaded  (open).  5/16  inch 
diameter  holes. 

The  above  examples  clearly  indicate  the  need  for  the  progressive  ply 
failure  procedure  ineorporated  in  the  SASCJ  code.  The  one-step  laminate 
level  failure  prediction  leads  to  the  dependence  of  the  failure  parameters  on 
the  laminate  layup.  Tlie  one-step  ply  level  failure  prediction  requires  the 
modification  of  actual  lamina  properties.  The  progressive  failure  procedure 
in  SASCJ  provides  a  meand  for  using  invariant  failure  parameters  and  ply 
property  degradation  rates  to  predict  the  strength  of  bolted  laminates.  An 
example  is  presented  below. 
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Figure  5-4.  Failure  Parameters  for  the  Laminate  Level  Prediction  of  Static 
Compressive  Strengths  of  Laminates  with  Unloaded  (Open)  5/16 
Inch  Diameter  Holes, 
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X, .  X^  and  S  (O’  ply  )  •  230,  267  and  t7.3  kal.  raapacllvaly 
X  k^and  S  ( i4S”  ply  }  >100,  100  and  100  kal.  raapdclkraly  (moddlad) 
X| ,  X^and  S  (OO"  ply  )  •  OS,  38.0  and  17.3  kal.  raapacllvaly  (modKiad) 


FlRure  5r-5.  Failure  Paraitieters  for  the  Ply  Level  Prediction  of  Static 
Tensile  Strengths  of  Laminates  with  Unloaded  (Open),  5/16 
Inch  Diameter  Holes, 
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Figure  3-6,  Failure  Parameters  for  the  Ply  Level  Prediction  of  Static 
Compressive  Strengths  of  Laminates  with  Unloaded  (Open), 
5/16  Inch  Diameter  Holes. 
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b.J  /i' t  ('TT!'. in  It  ion  of  i  >1  i  1  vir »'  rarnmot  rrr 

Coii^sldcr  a  lanundtt'  that  transfers  an  applied  tensile  or  compressive 

lead  in  double  shear  via  a  fastener  (see  Figure  5-7).  The  strength  of  this 

lai-'inate  is  predicted  using  the  progressive  ply  failure  proceduri'  in  SASCJ, 

assuning  bilinear  ply  behavior  (Figure  3-8)  and  average  stress  failure 

criteria.  The  strengths  of  the  +45^  plies  are  computed  using  lamination 

theory.  For  chch  ply  failure  mode  (net  section,  shear-out  and  bearing),  a 

characteristic  distance  (a  ,  a  or  a  .  ),  a  k^/k,  value  and  a 

initial  *21 

P  value  have :  to  be  determined.  Figure  5-7  presents  the  nine 

parameters  that  yielded  SASCJ  predictions  that  agreed  with  the  exper  ir ou la 1 
obsein/atlons  and  measurements  in  Reference  5-1,  including  laminate  strength 
and  failure  modes.  Note  that  only  one  failure  parameter  (a  )  had  to  be 
changed  to  make  predictions  correlate  accurately  with  measurements  corres¬ 
ponding  to  the  30/60/10  laminate.  If  the  a  value  for  the  tensile-loaded 

.  ons 

30/60/10  laminate  is  not  changed  from  the  0.10  inch  value,  only  a  error 
will  result  in  the  predicted  strength  (see  Section  5. A). 

Failure  parameters  based  on  a  linear  ply  behavior  (k2/kj  ■  0,  p'^itlnate^ 
P  “  1)  and  the  Tsai-Hill  criterion  or  the  maximum  strain  criterion 

arc  also  p.  asented  in  Figyre  5-7.  These  values,  corresponding  to  the  fully- 
loaded  hole  situation,  are  seen  to  be  different  from  those  corresponding 
to  the  unloaded  hole  situation  in  Figure  5-5.  They  are  also  very  dependent 
on  the  l.iralnate  layup. 

Based  on  the  above  results,  it  is  concluded  that  the  progressive  failure 
procedure  using  the, average  stress  failure  criterion  and  bilinear  ply  be¬ 
havior  is  the  best  choice  among  the  studied  methods.  Therefore,  subsequent 
correlation  studies  are  performed  using  this  procedure  with  invariant  failure 
parameters  and  ply  degradation  rates.  The  invariant  properties  are  obtained 
based  on  Imposed  correlation  between  SASCJ  predictions  and  test  measurements 
corresponding  to  a  double  shear  (full-bearing)  load  transfer  configuration. 

5* A  SASCJ  Strength  Predictions  for  Laminates  with  Fully-Loaded  Holes 
in  a  Single  Shear  Load  Transfer  Configuration 

For  ASl/3501-6  graphlte/opoxy  l«i>lnates,  tested  in  Reference  5-1,  the 
various  failure  parameters  obtained  through  a  trial  and  error  procedure  are 
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Figure  5-7.  Characteristic  Distances  and  Ply  Degradation  Parameters  for 
Laminates  with  Fully-Loaded  5/16  Inch  Diameter  Holes, 
Subjected  to  Static  Tensile  and  Compressive  Loading  In 
Double  Shear, 


0.10 


listed  in  Fifiure  5-7 .  These  include  the  characteristic  distances  (a  = 

ons 

inch,  a  ,  =  0.025  inch,  and  a  =0.08  inch),  the  ply  modulus  change 

corresponding  to  every  failure  mode  (k,/kj  =0.1  for  net  section,  bearing  and 
shear  out),  and  the  ultimate-to-inltial  ply  failure  load  ratio  for  each 
failure  mode  (1. 02  and  1 1 , 50  for  net  section  and  bearing,  and  1.12  and  1.80 
for  shear-out  under  tension  and  compression,  respectively).  The 
pinitial  ratio,  corresponding  to  the  shear-out  mode  of  ply  failure,  increases 
when  the  applied  load  is  cFianged  from  tension  to  compression,  due  to  the 
larger  shear-out  (or  shear-in)  area  under  this  loading  condition.  But,  for 
subsequent  correlation  studies,  is  assumed  to  be  1.12  for 

the  shear-out  mode  of  ply  failure,  irrespective  of  the  loading  mode.  This 
results  in  the  use  of  constant  failure  parameters  when  strength  predictions 
are  made  for  bolted  ASl/35Ql-6graphlte/epoxy  laminates  using  the  SASCJ  code. 

The  SASGJ  code  and  the  above  failure  parameters  (obtained  through  a 
forced  correlation  corresponding  to  a  double  shear  load  transfer  situation) 
were  subsequently  used  to  predict  the  strength  of  laminates  with  fully-loaded 
holes  in  a  single  shear  load  transfer  configuration.  Figures  5-8  and  5-9 
present  a  comparison  between  SASCJ  predictions  and  test  results  from 
Reference  5»Ii  With  the  exception  of  two  test  cases,  an  excellent  agreement 
between; SASCJ  predictions  and  test  results  is  noticed.  Under  tensile  loading, 
a  17%  erroi"  between  test  and  analysis  is  seen  for  30/60/10  (percentages  of 
0°,  +A5°  and  90®  plies,  respectively)  laminates  in  a  double  shear  configura¬ 
tion.  Under  compressive  loading,  a  27%  difference  between  the  SASCJ  pre- 
dictl''n  and  the  test  result  in  Reference  5-1,  for  70/20/10  laminates  in 
single  shear,  is  noticed.  Also,  a. 23%  difference  between  test  and  analysis 
is  seen  for  50/40/10  laminates  in  a  double  shear  load  transfer  configuration. 
These  differences,  under  compressive  loading,  will  be  reduced  if 
pinitial  shear-out  ply  failure  mode  is  assumed  to  vary  with  the  shear- 

out  area  (E/D).  This  will  be  attempted  later  when  the  SASCJ  code  is  used 
in  the  development  of  a  design  guide, 

5,5  Sample  SASCJ  Input/Output 

The  details  of  the  input  data  for  SASCJ  and  its  prediction  of  pro¬ 
gressive  ply  failures  are  presented  by  conalderlng  test  case  2  in  Figure  5-8. 


CM/I^  LAMIMATi  ' 


DOUBLE  SHEAR 


SINOLE  SHEAR 


LAVUR  r 

(*  0®.  146® 

TEST 

FAILURE  LOAD  (KIRS) 

TEST 

FAILURE  LOAD  (KIR 

AND  S0®RLIES)|  1 

(REF.  6-1) 

1  EXPT. 

(REF.  6-1) 

ANAL.*  ' 

CASE 

<ftEF.  9-1) 

1  EXPT. 
(R£F.  6-1) 

ANAL. 

60/40/10 

S2 

4.76 

4.67 

a 

4. SI 

4.66 

70/20/ 10 

62 

a. 9a 

1 

2.64 

a  1 

4.06 

4.04 

20/60/10 

64 

6. 27 

4.67 

aa 

•  .16 

6.0* 

*  *obrg  '  •••0  *  »S«  UrMp# 

E|/k,  ■  0.1  lor  •iiUllur*  med«« 

p  «S  ,  1.0S.  I.SO  SAd  1.1#  tdf  SfS  dAd  id,  r#»S»etl*dl» 


Figure  5-8,  SASCw^  Predictions  for  Sample  Laminates  with  Fully-Loaded 
5/16  Inch  Diameter  Holes,  Subjected  to  Static  Tensile 
Loading , 
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Figure  5-9^  SASCJ  Predictions  for  Sample  Laminates  with  Fully-Loaded 
S/16  Inch  pidmeter  Holes,  Subjected  to  Static  Compressive 
Loading, 
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In  this  sample  situation,  a  static  tensile  load  is  transferred  from  a  laminate 
to  aluminum  in  single  shear,  via  a  5/16  inch  diameter  steel  fastener.  The 
laminate  is  made  up  of  ASl/3501-6  graphlte/cpoxy,  and  has  a  50/40/10  layup 
defined  by  [(45/0/-45/0) 2  /0/90]g.  Its  E/D  and  W/D  values  are  3  and  6, 
respectively.  Table  5-2  presents  the  SASCJ  input  data  for  this  test  case. 

The  tensile  and  compressive  strengths,  and  the  inplane  shear  strengths,  for 
the  four  ply  types  (vd.th  fiber  orientations  of  0°,  45°,  -45°  and  90°)  were 
obtained  using  available  ASl/3501-6  graphlte/epoxy  properties  and  laminate 
theory.  Referring  to  Figure  3-8,  the  bilinear  load  versus  deflection  be¬ 
havior  was  defined  by  assuming  k2/kj  "0.1  for  all  the  failure  modes,  and 
pUltlmate^pinitial  assumed  to  be  1.02,  1.50  and  1.12  for  the  net  section, 

bearing  and  shear-out  modes  of  ply  failure,  respectively.  The  rotational 

constraints  at  the  fastener  head  and  nut  locations  were  represented  by 

12 

assuming  and  R2  values  to  be  10  (see  Equation  3-16). 

Table  5-3  presents  the  results  predicted  by  the  SASCJ  code  for  the 
problem  defined  above.  Nodes  1  to  30  are  within  the  aluminum  plate,  and 
nodes  31  to  70  are  within  the  laminate,  SASCJ  predicts  a  net  section  failure 
in  the  90°  piles  to  be  the  first  failure.  This  Is  followed  by  shear-out 
failures  In  -45°  piles,  net  section  failures  In  the  +45°  piles,  total  failures 
of  the  90°  piles,  shear-out  failures  In  the  0°  piles,  etc.  At  4.744  kips, 
total  failure  of  the  laminate  Is  predicted.  The  predicted  failure  load  Is 
within  5%  of  the  average  measured  value  of  4.914  kips  in  Reference  5-1,  Also, 
Table  5-3  Indicates  that  laminate  failure  is  predominately  precipitated  by 
shear-out  failure  in  the  0°  and  -45°  plies,  and  Is  accompanied  by  net  section 
failures  In  the  90°  and  +45°  plies.  The  prlaiary  mode  of  failure  observed  in 
Reference  5-1  Is  the  predicted  "partial  shear-out"  of  the  laminate. 

5 . 6  SASCJ  Strength  Predictions  for  Sample  Laminates  with  Partially- 

Loaded  Holes  —  The  General  Single  Fastener  Situation 

Laminates  with  partially-loaded  holes  are  analyzed  by  superimposing 
solutions  for  the  corresponding  unloaded  (open)  and  fully-loaded  hole  prob¬ 
lems  (see  Section  4.3),  Section  5.3  discusses  how  the  various  failure  para¬ 
meters  are  obtained.  Using  these  parameters,  SASCJ  Is  used  to  predict  the 
strength  of  three  laminates  subjected  to  partial  fastener  loads.  Figure  5-10 
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Figure  5-10.  Comparison  of  SASCJ  Predictions  with  Experimental 
Measurements  for  50/40/10  ASl/3501-6  Laminates 
vrlth  Partially-Loaded  5/16  Inch  Diameter  Holes, 
Under  Tensile  Loading. 
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compares  SASCJ  predictions  with  experimental  measurements  (Reference  5-1) 

%^en  a  tensile  load  is  applied  on  50/40/10,  ASl/3501-6,  graphlte/epoxy  lami¬ 
nates  with  partially-loaded  5/16  inch  diameter  holes.  SASCJ  predictions 
agree  very  well  with  experimental  measurements  for  the  listed  fastener 
bearing-to-total  load  ratios. 

Similar  results  for  tension-loaded  70/20/10  laminates,  with  partially- 
loaded  5/16  inch  diameter  holes,  are  presented  In  Figure  5-11.  The  maximum 
difference  between  SASCJ  predictions  and  test  results  for  this  case  is  seen 
to  be  approximately  30%.  A  better  correlation,  between  test  and  analysis  is 
demonstrated  by  tension-loaded  30/60/10  lamina tea  with  partially-loaded  holes 
fsee  Figure  5-12). 

Figures  5-13  to  5-15  compare  SASCJ  predictions  with  test  results  from 
Reference  5-1,  corresponding  to  a  compressive  applied  load  on  50/40/10, 
70/20/10.  and  30/60/10  laminates  with  partially  loaded,  5/16  inch  diameter 
holes.  Again,  SASCJ  predictions  agree  very 'Well  with  test  results  on 
50/40/10  and  30/60/10  laminates,  and  agree  reasonabty  well  with  test  results 
on  70/20/10  laminates,  fori  the  listed  bearing/total  load  ratios. 

Figures  5-10  to  S-lStdemonstrate  an  acceptable  correlation  between 
analytical  predictions  and  test  results  for  a  general  single  fastener  loca¬ 
tion  in  a  bolted  laminate.  This,  In  conjunction  with  the  demonstrated 
correlation  In  Section  5.4,  validates  the  use  of  SASCJ  for  predicting  the 
strength  of  bolted  laminates. 
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Figure  5-11.  Comparison  of  SASCJ  Predictions  with  Experimental 
Measurements  for  70/20/10  ASl/3501-6  Laminates 
with  Partially-Loaded  5/16  Inch  Diameter  Holes, 
Under  Tensile  Loading, 
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Figure  5-.13,  Comparison  of  SASCJ  Predictions  with  Experimental 

Measurements  for  50/40/10  ASl/ 3501-6  Laminates  with 
Partially-Loaded  5/16  Inch  Diameter  Holes,  Under 
Compressive  Loading, 
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Figure  5-15,  Comparison  of  SASCJ  Fred lotions with  Experimental 
Measurements  for  30/60/10  ASl/3501-6  Laminates 
with  Partially-Loaded  5/16  Inch  Diameter  Holes, 
under  Compressive  Loading. 
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'  SECTION  6  ' 

CONCLUSIONS 

(SASCJ)  was  developed  for  laminated  and/or  metallic 
plates  bolted  together  by  a  single  fastener.  This  analysis  includes  a  two-  ^ 

dimensional  analysis  of  a  finite  anisotropic  plate  with  a  loaded  or  unloaded 
hole  (FIGECM),  a  fastener  analysis  that  enables  an  approximate  prediction  of 
the  local  three  dimensional  stress  state,  accounting  for  fastener  flexibility, 
load  eccentricity,  fastener  torque,  etc.  A),  and  a  progressive  failure 
procedure  based  on  bilinear  ply  behavior.  FIGECM  and  FDFA  predictions  were 
validated  using  available  numerical  and  analytical  solutions.  SASCJ  strength 
predictions  were  shown  to  correlate  very  well  with  test  results  from 
Reference  5-1. 

The  major  accomplishments  of  the  reported  program  cask  are  sunmarlzed 

(1)  The  developed  SASCJ  code  accounts  for  finite  (actual)  planform  dlmen- 

of  bolted  plates.  It  is,  therefore,  capable  of  predicting  the 
:  E/D  and  W/D  values  (<  3  and  6,  respectlve3.y)  on  the 

strength  of  bolted  laminates.  , 

(2)  The  combination  of  the  two-dimensional  anisotropic  plate  analysis  ^ 

^  approximately  accounts  for 

the  complex  three-dimensional  stress/strain  state  at  the  fastener  loca¬ 
tion.  This  also  enables  an  analytical  assessment  of  the  effect  of  im¬ 
portant  joint  parameters  that  Include  fastener  flexibility,  load  / 

eccentricity  (single  versus  double  shear  load  transfer),  fastener  torque, 
laminate  stacking  sequence,  etc. 

(3)  The  progressive  failure  procedure  in  the  SASCJ  code,  incorporating  a  bi¬ 
linear  ply  behavior  into  average  stress  failure  criteria,  provides  sig¬ 
nificant  improvements  in  the  strength  predictive  capability.  It  en¬ 
ables  the  use  of  unmodified  ply  properties  and  invariant  failure  para¬ 
meters  to  predict  the  sequence  of  ply  failures  until  one  of  the  bolted 
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plates  suffers  total  failure  (joint  failure).  It  predicts  the  failure 
modes  at  the  ply  level,  accounting  for  a  realistic  combination  of  many 
failure  modes  in  the  bolted  laminate. 

(4)  The  invariant  failure  parameters,  required  as  input  by  the  SASCJ  code, 
are  obtainable  from  a  limited  number  of  tests  on  bolted  laminates  in  a 
double  shear  load  transfer  configuration.  These  include  the  character¬ 
istic  distances  for  the  average  stress  criteria  and  the  k2/kj  and 

_ultlmate/_lnltlal  r  ji. 

P  P  ratios  for  net  section,  bearing  and  shear-out 

failure  modes.  This  report  Includes  failure  parameter  values  for  the 
ASl/3501-6  graphite/epoxy  material  system.  For  a  different  material 
system,  only  a  few  double  shear,  full-bearing  tests  have  to  be  con¬ 
ducted  to  determine  the  applicable  failure  parameters. 

(5)  The  SASCJ  code  can  predict  the  strength  of  an  Isolated  single  fastener 
region  in  a  laminate  bolted  to  another  plate  using  many  fasteners,  pro¬ 
vided  the  fastener  load  is  known.  The  SASCJ  code  can,  therefore,  be 
directly  combined  with  a  two-dimensional  load  distribution  analysis 

to  predict  the  strength  of  laminates  bolted  to  other  structural  elements 
using  many  fasteners.  This  strength  analysis  of  multifastener  joints  In 
composite  structures  Is  currently  under  development  In  the  ongoing 
Northrop/AFWAL  program. 

(6)  The  SASCJ  code  Is  currently  being  used  to  develop  a  design  guide  for 
bolted  joints.  Expressions  are  being  developed  to  quantify  the  effect 
of  fastener  torque  on  the  joint  strength.  The  effects  of  the  geometry 
and  the  properties  of  the  bolted  plates  and  the  fastener  are  being 
systematically  studied.  Obtained  results  are  being  analyzed  and  trans¬ 
lated  Into  applicable  design  guidelines. 

The  above  accomplishments  of  the  SASCJ  code  represent  a  significant 
Improvement  over  the  currently  available  strength  analysis  code  (BJSFM). 
Nevertheless,  SASCJ  does  contain  a  restriction  and  a  limitation  that  should 
be  addressed  in  future  developmental  efforts.  The  code  Is  currently  re¬ 
stricted  to  a  protruding  head  fastener  geometry,  and  cannot  predict  the 
effect  of  a  flush-head  (countersunk)  fastener  geometry.  Also,  the  fastener 
bearing  load  Is  always  assumed  to  be  cosinusoidal  in  form,  and  Is  assumed 
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to  act  over  half  the  hole  boundary.  This  results  in  a  state  of  stress  that 

near  the  hole  boundary.  The  computed  ply  stresses  are, 
therefore,  accurate  only  at  planform  locations  that  are  at  least  a  laminate 
thickness  away  from  the  hole  boundary.  But,  failure  predictions  are  based  on 
stresses  computed  oyer  characteristic  distances  that  are  smaller  or  nearly 
equal  in  magnitude  compared  to  the  laminate  thickness.  This  leads  to  an 
inaccuracy  in  the  predicted  joint  strength,  and  the  degree  of  inaccuracy  is 
unknown.  A  modification  of  the  FIGEOM  and  FDFA  analyses  in  tie  SASCJ  code 
is  strongly  recommended  as  a  future  analytical  effort  to  eliminate  the 
mentioned  restriction  and  limitation. 
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